martinEcarnot/signe
Useful functions to handle spectra

Package index
Search the martinEcarnot/signe package
Vignettes
Functions
0
Source code
78
Man pages
0
Home
/
GitHub
/
martinEcarnot/signe
/
Script_R_2020/Non utilise/Doublons d'un programmes actuel/PLSDA_souches_cepages.R

Script_R_2020/Non utilise/Doublons d'un programmes actuel/PLSDA_souches_cepages.R
In martinEcarnot/signe: Useful functions to handle spectra 

library(MASS)
#library(mixOmics)
library(FactoMineR)
library(signal)
library(plyr)
library(caret)
library(dplyr)
library(rnirs)

rm(list = ls())

source('C:/Users/avitvale/Documents/Script_R/adj_asd.R')
source('C:/Users/avitvale/Documents/Script_R/SIGNE_load.R')
source('C:/Users/avitvale/Documents/Script_R/SIGNE_maha.R')
source('C:/Users/avitvale/Documents/Script_R/SIGNE_maha0.R')


# Choix de la fixation du tirage aleatoire (pour comparaison, rend les repetitions inutiles)
#set.seed(1)
brb4="C:/Users/avitvale/Documents/Test/globalmatrix"
load(file=brb4)

## Filtrage des spectres aberrants
globalmatrix=globalmatrix[globalmatrix[,500]>0.6,]
globalmatrix=globalmatrix[globalmatrix[,1]<0.2,]
globalmatrix=globalmatrix[globalmatrix[,2000]<0.25,]

# Data Filter
# Select dates
dates=list(
#   "20180619T"
# "20180619P"
 # "20180619N"
 #  ,"20180627T"
# , "20180627P"
#   ,"20180627N"
 #"20180704P"
  # "20180709T"
 # "20180709P",
#  "20180709N"
#   ,"20180710P"
  # ,"20180816T"
  #, "20180816P"
  # ,"20180816N"
   #"20180817T"
  # ,"20180817P"
 # ,"20180817N"
  # "20170524P"
  # , "20170529P"
   # ,"20170606P"
 #   ,"20170612P"
   "20180619P"
 #    ,"20170626P"
 # ,  "20170703P"
 # ,"20170710P"
 #  ,"20170717P"
 #  ,"20170724P"
 #  ,"20170731P"
# ,"20180823P"
  # ,"20180823T"
  # ,"20180823N"
 #  ,"20180731P"
  # ,"20180731T"
  # ,"20180731N"
 #  ,"20180810P"
  # ,"20180810T"
  # ,"20180810N"
 #  ,"20180724P"
  # ,"20180724T"
  # ,"20180724N"
 # "20180809T"
)

iok=substr(rownames(globalmatrix),1,9) %in% dates
globalmatrix=globalmatrix[iok,]

### FIXATION DES PARAMETRES UTILISES:
## Nombre de repetitions de la boucle de PLSDA:
repet= 50
## Parametres du Savitsky-Golay (p=degre du polynome, n= taille de la fenetre, m=ordre de derivation)
p=2
n=11
m=1
## Nombre de VL max autorisees
ncmax=10
## Taille de l'echantillon de validation (1/v):
#v=3
## Nombre de groupes de CV
k=2

## PLSDA ##
sp=globalmatrix

### Pretraitements
## Ajustement des sauts de detecteurs (Montpellier: sauts ?? 1000 (651 eme l.o.) et 1800 (1451))
sp=adj_asd(sp,c(602,1402))
## Reduction des variables (extremites bruitees)
# sp=sp[,seq(51,ncol(sp)-30,1)]
## SNV
sp=t(scale(t(sp)))
## Derivation Savitsky Golay
sp=t(apply(sp,1,sgolayfilt,p=p,n=n,m=m))



iok2=substr(rownames(sp),1,9) %in% dates
sp=sp[iok2,]

## Creation de la matrice de classes
class=as.factor(substr(rownames(sp),9,9))
## Variable qui mesure le nombre de classes
c=length(levels(class))

###Echantillonnage par "souches"
##Separation de sp en 30 matrices
titre=rownames(sp)

# Boucle en fonction de la date
# d1 = (substr(titre,1,9))== "20170624P"
# d2 = (substr(titre,1,9))== ""
#
# #g1=sp[(d1==TRUE),] #1
#
#
# s=2 #indice de la date
#
# for (i in 1:s) {
#   commd=paste("g",i,"=sp[(d",i,sep="")
#   eval(parse(text=commd))
#
# }

#Tous
x1= (substr(titre,11,16))== "015-01" |  (substr(titre,11,16))== "015-02" |  (substr(titre,11,16))== "015-03" |  (substr(titre,11,16))== "015-04" |  (substr(titre,11,16))== "015-05" |  (substr(titre,11,16))== "015-06"
x11=(substr(titre,11,16))== "015-07" |  (substr(titre,11,16))== "015-08" |  (substr(titre,11,16))== "015-09" |  (substr(titre,11,16))== "015-10" |  (substr(titre,11,16))== "015-11" |  (substr(titre,11,16))== "015-12"
x21=(substr(titre,11,16))== "015-13" |  (substr(titre,11,16))== "015-14" |  (substr(titre,11,16))== "015-15" |  (substr(titre,11,16))== "015-16" |  (substr(titre,11,16))== "015-17" |  (substr(titre,11,16))== "015-18"
x2= (substr(titre,11,16))== "169-01" |  (substr(titre,11,16))== "169-02" |  (substr(titre,11,16))== "169-03" |  (substr(titre,11,16))== "169-04" |  (substr(titre,11,16))== "169-05" |  (substr(titre,11,16))== "169-06"
x12=(substr(titre,11,16))== "169-07" |  (substr(titre,11,16))== "169-08" |  (substr(titre,11,16))== "169-09" |  (substr(titre,11,16))== "169-10" |  (substr(titre,11,16))== "169-11" |  (substr(titre,11,16))== "169-12"
x22=(substr(titre,11,16))== "169-13" |  (substr(titre,11,16))== "169-14" |  (substr(titre,11,16))== "169-15" |  (substr(titre,11,16))== "169-16" |  (substr(titre,11,16))== "169-17" |  (substr(titre,11,16))== "169-18"
x3= (substr(titre,11,16))== "222-01" |  (substr(titre,11,16))== "222-02" |  (substr(titre,11,16))== "222-03" |  (substr(titre,11,16))== "222-04" |  (substr(titre,11,16))== "222-05" |  (substr(titre,11,16))== "222-06"
x13=(substr(titre,11,16))== "222-07" |  (substr(titre,11,16))== "222-08" |  (substr(titre,11,16))== "222-09" |  (substr(titre,11,16))== "222-10" |  (substr(titre,11,16))== "222-11" |  (substr(titre,11,16))== "222-12"
x23=(substr(titre,11,16))== "222-13" |  (substr(titre,11,16))== "222-14" |  (substr(titre,11,16))== "222-15" |  (substr(titre,11,16))== "222-16" |  (substr(titre,11,16))== "222-17" |  (substr(titre,11,16))== "222-18"
x4= (substr(titre,11,16))== "471-01" |  (substr(titre,11,16))== "471-02" |  (substr(titre,11,16))== "471-03" |  (substr(titre,11,16))== "471-04" |  (substr(titre,11,16))== "471-05" |  (substr(titre,11,16))== "471-06"
x14=(substr(titre,11,16))== "471-07" |  (substr(titre,11,16))== "471-08" |  (substr(titre,11,16))== "471-09" |  (substr(titre,11,16))== "471-10" |  (substr(titre,11,16))== "471-11" |  (substr(titre,11,16))== "471-12"
x24=(substr(titre,11,16))== "471-13" |  (substr(titre,11,16))== "471-14" |  (substr(titre,11,16))== "471-15" |  (substr(titre,11,16))== "471-16" |  (substr(titre,11,16))== "471-17" |  (substr(titre,11,16))== "471-18"
x5= (substr(titre,11,16))== "509-01" |  (substr(titre,11,16))== "509-02" |  (substr(titre,11,16))== "509-03" |  (substr(titre,11,16))== "509-04" |  (substr(titre,11,16))== "509-05" |  (substr(titre,11,16))== "509-06"
x15=(substr(titre,11,16))== "509-07" |  (substr(titre,11,16))== "509-08" |  (substr(titre,11,16))== "509-09" |  (substr(titre,11,16))== "509-10" |  (substr(titre,11,16))== "509-11" |  (substr(titre,11,16))== "509-12"
x25=(substr(titre,11,16))== "509-13" |  (substr(titre,11,16))== "509-14" |  (substr(titre,11,16))== "509-15" |  (substr(titre,11,16))== "509-16" |  (substr(titre,11,16))== "509-17" |  (substr(titre,11,16))== "509-18"
x6= (substr(titre,11,16))== "525-01" |  (substr(titre,11,16))== "525-02" |  (substr(titre,11,16))== "525-03" |  (substr(titre,11,16))== "525-04" |  (substr(titre,11,16))== "525-05" |  (substr(titre,11,16))== "525-06"
x16=(substr(titre,11,16))== "525-07" |  (substr(titre,11,16))== "525-08" |  (substr(titre,11,16))== "525-09" |  (substr(titre,11,16))== "525-10" |  (substr(titre,11,16))== "525-11" |  (substr(titre,11,16))== "525-12"
x26=(substr(titre,11,16))== "525-13" |  (substr(titre,11,16))== "525-14" |  (substr(titre,11,16))== "525-15" |  (substr(titre,11,16))== "525-16" |  (substr(titre,11,16))== "525-17" |  (substr(titre,11,16))== "525-18"
x7= (substr(titre,11,16))== "685-01" |  (substr(titre,11,16))== "685-02" |  (substr(titre,11,16))== "685-03" |  (substr(titre,11,16))== "685-04" |  (substr(titre,11,16))== "685-05" |  (substr(titre,11,16))== "685-06"
x17=(substr(titre,11,16))== "685-07" |  (substr(titre,11,16))== "685-08" |  (substr(titre,11,16))== "685-09" |  (substr(titre,11,16))== "685-10" |  (substr(titre,11,16))== "685-11" |  (substr(titre,11,16))== "685-12"
x27=(substr(titre,11,16))== "685-13" |  (substr(titre,11,16))== "685-14" |  (substr(titre,11,16))== "685-15" |  (substr(titre,11,16))== "685-16" |  (substr(titre,11,16))== "685-17" |  (substr(titre,11,16))== "685-18"
x8= (substr(titre,11,16))== "747-01" |  (substr(titre,11,16))== "747-02" |  (substr(titre,11,16))== "747-03" |  (substr(titre,11,16))== "747-04" |  (substr(titre,11,16))== "747-05" |  (substr(titre,11,16))== "747-06"
x18=(substr(titre,11,16))== "747-07" |  (substr(titre,11,16))== "747-08" |  (substr(titre,11,16))== "747-09" |  (substr(titre,11,16))== "747-10" |  (substr(titre,11,16))== "747-11" |  (substr(titre,11,16))== "747-12"
x28=(substr(titre,11,16))== "747-13" |  (substr(titre,11,16))== "747-14" |  (substr(titre,11,16))== "747-15" |  (substr(titre,11,16))== "747-16" |  (substr(titre,11,16))== "747-17" |  (substr(titre,11,16))== "747-18"
x9= (substr(titre,11,16))== "787-01" |  (substr(titre,11,16))== "787-02" |  (substr(titre,11,16))== "787-03" |  (substr(titre,11,16))== "787-04" |  (substr(titre,11,16))== "787-05" |  (substr(titre,11,16))== "787-06"
x19=(substr(titre,11,16))== "787-07" |  (substr(titre,11,16))== "787-08" |  (substr(titre,11,16))== "787-09" |  (substr(titre,11,16))== "787-10" |  (substr(titre,11,16))== "787-11" |  (substr(titre,11,16))== "787-12"
x29=(substr(titre,11,16))== "787-13" |  (substr(titre,11,16))== "787-14" |  (substr(titre,11,16))== "787-15" |  (substr(titre,11,16))== "787-16" |  (substr(titre,11,16))== "787-17" |  (substr(titre,11,16))== "787-18"
x10=(substr(titre,11,16))== "877-01" |  (substr(titre,11,16))== "877-02" |  (substr(titre,11,16))== "877-03" |  (substr(titre,11,16))== "877-04" |  (substr(titre,11,16))== "877-05" |  (substr(titre,11,16))== "877-06"
x20=(substr(titre,11,16))== "877-07" |  (substr(titre,11,16))== "877-08" |  (substr(titre,11,16))== "877-09" |  (substr(titre,11,16))== "877-10" |  (substr(titre,11,16))== "877-11" |  (substr(titre,11,16))== "877-12"
x30=(substr(titre,11,16))== "877-13" |  (substr(titre,11,16))== "877-14" |  (substr(titre,11,16))== "877-15" |  (substr(titre,11,16))== "877-16" |  (substr(titre,11,16))== "877-17" |  (substr(titre,11,16))== "877-18"


# #sans 747
# x1= (substr(titre,11,16))== "015-01" |  (substr(titre,11,16))== "015-02" |  (substr(titre,11,16))== "015-03" |  (substr(titre,11,16))== "015-04" |  (substr(titre,11,16))== "015-05" |  (substr(titre,11,16))== "015-06"
# x10=(substr(titre,11,16))== "015-07" |  (substr(titre,11,16))== "015-08" |  (substr(titre,11,16))== "015-09" |  (substr(titre,11,16))== "015-10" |  (substr(titre,11,16))== "015-11" |  (substr(titre,11,16))== "015-12"
# x19=(substr(titre,11,16))== "015-13" |  (substr(titre,11,16))== "015-14" |  (substr(titre,11,16))== "015-15" |  (substr(titre,11,16))== "015-16" |  (substr(titre,11,16))== "015-17" |  (substr(titre,11,16))== "015-18"
# x2= (substr(titre,11,16))== "169-01" |  (substr(titre,11,16))== "169-02" |  (substr(titre,11,16))== "169-03" |  (substr(titre,11,16))== "169-04" |  (substr(titre,11,16))== "169-05" |  (substr(titre,11,16))== "169-06"
# x11=(substr(titre,11,16))== "169-07" |  (substr(titre,11,16))== "169-08" |  (substr(titre,11,16))== "169-09" |  (substr(titre,11,16))== "169-10" |  (substr(titre,11,16))== "169-11" |  (substr(titre,11,16))== "169-12"
# x20=(substr(titre,11,16))== "169-13" |  (substr(titre,11,16))== "169-14" |  (substr(titre,11,16))== "169-15" |  (substr(titre,11,16))== "169-16" |  (substr(titre,11,16))== "169-17" |  (substr(titre,11,16))== "169-18"
# x3= (substr(titre,11,16))== "222-01" |  (substr(titre,11,16))== "222-02" |  (substr(titre,11,16))== "222-03" |  (substr(titre,11,16))== "222-04" |  (substr(titre,11,16))== "222-05" |  (substr(titre,11,16))== "222-06"
# x12=(substr(titre,11,16))== "222-07" |  (substr(titre,11,16))== "222-08" |  (substr(titre,11,16))== "222-09" |  (substr(titre,11,16))== "222-10" |  (substr(titre,11,16))== "222-11" |  (substr(titre,11,16))== "222-12"
# x21=(substr(titre,11,16))== "222-13" |  (substr(titre,11,16))== "222-14" |  (substr(titre,11,16))== "222-15" |  (substr(titre,11,16))== "222-16" |  (substr(titre,11,16))== "222-17" |  (substr(titre,11,16))== "222-18"
# x4= (substr(titre,11,16))== "471-01" |  (substr(titre,11,16))== "471-02" |  (substr(titre,11,16))== "471-03" |  (substr(titre,11,16))== "471-04" |  (substr(titre,11,16))== "471-05" |  (substr(titre,11,16))== "471-06"
# x13=(substr(titre,11,16))== "471-07" |  (substr(titre,11,16))== "471-08" |  (substr(titre,11,16))== "471-09" |  (substr(titre,11,16))== "471-10" |  (substr(titre,11,16))== "471-11" |  (substr(titre,11,16))== "471-12"
# x22=(substr(titre,11,16))== "471-13" |  (substr(titre,11,16))== "471-14" |  (substr(titre,11,16))== "471-15" |  (substr(titre,11,16))== "471-16" |  (substr(titre,11,16))== "471-17" |  (substr(titre,11,16))== "471-18"
# x5= (substr(titre,11,16))== "509-01" |  (substr(titre,11,16))== "509-02" |  (substr(titre,11,16))== "509-03" |  (substr(titre,11,16))== "509-04" |  (substr(titre,11,16))== "509-05" |  (substr(titre,11,16))== "509-06"
# x14=(substr(titre,11,16))== "509-07" |  (substr(titre,11,16))== "509-08" |  (substr(titre,11,16))== "509-09" |  (substr(titre,11,16))== "509-10" |  (substr(titre,11,16))== "509-11" |  (substr(titre,11,16))== "509-12"
# x23=(substr(titre,11,16))== "509-13" |  (substr(titre,11,16))== "509-14" |  (substr(titre,11,16))== "509-15" |  (substr(titre,11,16))== "509-16" |  (substr(titre,11,16))== "509-17" |  (substr(titre,11,16))== "509-18"
# x6= (substr(titre,11,16))== "525-01" |  (substr(titre,11,16))== "525-02" |  (substr(titre,11,16))== "525-03" |  (substr(titre,11,16))== "525-04" |  (substr(titre,11,16))== "525-05" |  (substr(titre,11,16))== "525-06"
# x15=(substr(titre,11,16))== "525-07" |  (substr(titre,11,16))== "525-08" |  (substr(titre,11,16))== "525-09" |  (substr(titre,11,16))== "525-10" |  (substr(titre,11,16))== "525-11" |  (substr(titre,11,16))== "525-12"
# x24=(substr(titre,11,16))== "525-13" |  (substr(titre,11,16))== "525-14" |  (substr(titre,11,16))== "525-15" |  (substr(titre,11,16))== "525-16" |  (substr(titre,11,16))== "525-17" |  (substr(titre,11,16))== "525-18"
# x7= (substr(titre,11,16))== "685-01" |  (substr(titre,11,16))== "685-02" |  (substr(titre,11,16))== "685-03" |  (substr(titre,11,16))== "685-04" |  (substr(titre,11,16))== "685-05" |  (substr(titre,11,16))== "685-06"
# x16=(substr(titre,11,16))== "685-07" |  (substr(titre,11,16))== "685-08" |  (substr(titre,11,16))== "685-09" |  (substr(titre,11,16))== "685-10" |  (substr(titre,11,16))== "685-11" |  (substr(titre,11,16))== "685-12"
# x25=(substr(titre,11,16))== "685-13" |  (substr(titre,11,16))== "685-14" |  (substr(titre,11,16))== "685-15" |  (substr(titre,11,16))== "685-16" |  (substr(titre,11,16))== "685-17" |  (substr(titre,11,16))== "685-18"
# x8= (substr(titre,11,16))== "787-01" |  (substr(titre,11,16))== "787-02" |  (substr(titre,11,16))== "787-03" |  (substr(titre,11,16))== "787-04" |  (substr(titre,11,16))== "787-05" |  (substr(titre,11,16))== "787-06"
# x17=(substr(titre,11,16))== "787-07" |  (substr(titre,11,16))== "787-08" |  (substr(titre,11,16))== "787-09" |  (substr(titre,11,16))== "787-10" |  (substr(titre,11,16))== "787-11" |  (substr(titre,11,16))== "787-12"
# x26=(substr(titre,11,16))== "787-13" |  (substr(titre,11,16))== "787-14" |  (substr(titre,11,16))== "787-15" |  (substr(titre,11,16))== "787-16" |  (substr(titre,11,16))== "787-17" |  (substr(titre,11,16))== "787-18"
# x9=(substr(titre,11,16))== "877-01" |  (substr(titre,11,16))== "877-02" |  (substr(titre,11,16))== "877-03" |  (substr(titre,11,16))== "877-04" |  (substr(titre,11,16))== "877-05" |  (substr(titre,11,16))== "877-06"
# x18=(substr(titre,11,16))== "877-07" |  (substr(titre,11,16))== "877-08" |  (substr(titre,11,16))== "877-09" |  (substr(titre,11,16))== "877-10" |  (substr(titre,11,16))== "877-11" |  (substr(titre,11,16))== "877-12"
# x27=(substr(titre,11,16))== "877-13" |  (substr(titre,11,16))== "877-14" |  (substr(titre,11,16))== "877-15" |  (substr(titre,11,16))== "877-16" |  (substr(titre,11,16))== "877-17" |  (substr(titre,11,16))== "877-18"



# #Comparaison 747 et 877
# x1= (substr(titre,11,16))== "015-01" |  (substr(titre,11,16))== "015-02" |  (substr(titre,11,16))== "015-03" |  (substr(titre,11,16))== "015-04" |  (substr(titre,11,16))== "015-05" |  (substr(titre,11,16))== "015-06"
# x9= (substr(titre,11,16))== "015-07" |  (substr(titre,11,16))== "015-08" |  (substr(titre,11,16))== "015-09" |  (substr(titre,11,16))== "015-10" |  (substr(titre,11,16))== "015-11" |  (substr(titre,11,16))== "015-12"
# x17=(substr(titre,11,16))== "015-13" |  (substr(titre,11,16))== "015-14" |  (substr(titre,11,16))== "015-15" |  (substr(titre,11,16))== "015-16" |  (substr(titre,11,16))== "015-17" |  (substr(titre,11,16))== "015-18"
# x2= (substr(titre,11,16))== "169-01" |  (substr(titre,11,16))== "169-02" |  (substr(titre,11,16))== "169-03" |  (substr(titre,11,16))== "169-04" |  (substr(titre,11,16))== "169-05" |  (substr(titre,11,16))== "169-06"
# x10=(substr(titre,11,16))== "169-07" |  (substr(titre,11,16))== "169-08" |  (substr(titre,11,16))== "169-09" |  (substr(titre,11,16))== "169-10" |  (substr(titre,11,16))== "169-11" |  (substr(titre,11,16))== "169-12"
# x18=(substr(titre,11,16))== "169-13" |  (substr(titre,11,16))== "169-14" |  (substr(titre,11,16))== "169-15" |  (substr(titre,11,16))== "169-16" |  (substr(titre,11,16))== "169-17" |  (substr(titre,11,16))== "169-18"
# x3= (substr(titre,11,16))== "222-01" |  (substr(titre,11,16))== "222-02" |  (substr(titre,11,16))== "222-03" |  (substr(titre,11,16))== "222-04" |  (substr(titre,11,16))== "222-05" |  (substr(titre,11,16))== "222-06"
# x11=(substr(titre,11,16))== "222-07" |  (substr(titre,11,16))== "222-08" |  (substr(titre,11,16))== "222-09" |  (substr(titre,11,16))== "222-10" |  (substr(titre,11,16))== "222-11" |  (substr(titre,11,16))== "222-12"
# x19=(substr(titre,11,16))== "222-13" |  (substr(titre,11,16))== "222-14" |  (substr(titre,11,16))== "222-15" |  (substr(titre,11,16))== "222-16" |  (substr(titre,11,16))== "222-17" |  (substr(titre,11,16))== "222-18"
# x4= (substr(titre,11,16))== "509-01" |  (substr(titre,11,16))== "509-02" |  (substr(titre,11,16))== "509-03" |  (substr(titre,11,16))== "509-04" |  (substr(titre,11,16))== "509-05" |  (substr(titre,11,16))== "509-06"
# x12=(substr(titre,11,16))== "509-07" |  (substr(titre,11,16))== "509-08" |  (substr(titre,11,16))== "509-09" |  (substr(titre,11,16))== "509-10" |  (substr(titre,11,16))== "509-11" |  (substr(titre,11,16))== "509-12"
# x20=(substr(titre,11,16))== "509-13" |  (substr(titre,11,16))== "509-14" |  (substr(titre,11,16))== "509-15" |  (substr(titre,11,16))== "509-16" |  (substr(titre,11,16))== "509-17" |  (substr(titre,11,16))== "509-18"
# x5= (substr(titre,11,16))== "685-01" |  (substr(titre,11,16))== "685-02" |  (substr(titre,11,16))== "685-03" |  (substr(titre,11,16))== "685-04" |  (substr(titre,11,16))== "685-05" |  (substr(titre,11,16))== "685-06"
# x13=(substr(titre,11,16))== "685-07" |  (substr(titre,11,16))== "685-08" |  (substr(titre,11,16))== "685-09" |  (substr(titre,11,16))== "685-10" |  (substr(titre,11,16))== "685-11" |  (substr(titre,11,16))== "685-12"
# x21=(substr(titre,11,16))== "685-13" |  (substr(titre,11,16))== "685-14" |  (substr(titre,11,16))== "685-15" |  (substr(titre,11,16))== "685-16" |  (substr(titre,11,16))== "685-17" |  (substr(titre,11,16))== "685-18"
# x6= (substr(titre,11,16))== "747-01" |  (substr(titre,11,16))== "747-02" |  (substr(titre,11,16))== "747-03" |  (substr(titre,11,16))== "747-04" |  (substr(titre,11,16))== "747-05" |  (substr(titre,11,16))== "747-06"
# x14=(substr(titre,11,16))== "747-07" |  (substr(titre,11,16))== "747-08" |  (substr(titre,11,16))== "747-09" |  (substr(titre,11,16))== "747-10" |  (substr(titre,11,16))== "747-11" |  (substr(titre,11,16))== "747-12"
# x22=(substr(titre,11,16))== "747-13" |  (substr(titre,11,16))== "747-14" |  (substr(titre,11,16))== "747-15" |  (substr(titre,11,16))== "747-16" |  (substr(titre,11,16))== "747-17" |  (substr(titre,11,16))== "747-18"
# x7= (substr(titre,11,16))== "787-01" |  (substr(titre,11,16))== "787-02" |  (substr(titre,11,16))== "787-03" |  (substr(titre,11,16))== "787-04" |  (substr(titre,11,16))== "787-05" |  (substr(titre,11,16))== "787-06"
# x15=(substr(titre,11,16))== "787-07" |  (substr(titre,11,16))== "787-08" |  (substr(titre,11,16))== "787-09" |  (substr(titre,11,16))== "787-10" |  (substr(titre,11,16))== "787-11" |  (substr(titre,11,16))== "787-12"
# x23=(substr(titre,11,16))== "787-13" |  (substr(titre,11,16))== "787-14" |  (substr(titre,11,16))== "787-15" |  (substr(titre,11,16))== "787-16" |  (substr(titre,11,16))== "787-17" |  (substr(titre,11,16))== "787-18"
# x8= (substr(titre,11,16))== "877-01" |  (substr(titre,11,16))== "877-02" |  (substr(titre,11,16))== "877-03" |  (substr(titre,11,16))== "877-04" |  (substr(titre,11,16))== "877-05" |  (substr(titre,11,16))== "877-06"
# x16=(substr(titre,11,16))== "877-07" |  (substr(titre,11,16))== "877-08" |  (substr(titre,11,16))== "877-09" |  (substr(titre,11,16))== "877-10" |  (substr(titre,11,16))== "877-11" |  (substr(titre,11,16))== "877-12"
# x24=(substr(titre,11,16))== "877-13" |  (substr(titre,11,16))== "877-14" |  (substr(titre,11,16))== "877-15" |  (substr(titre,11,16))== "877-16" |  (substr(titre,11,16))== "877-17" |  (substr(titre,11,16))== "877-18"

# #Comparaison 471 et 525
# x1= (substr(titre,11,16))== "015-01" |  (substr(titre,11,16))== "015-02" |  (substr(titre,11,16))== "015-03" |  (substr(titre,11,16))== "015-04" |  (substr(titre,11,16))== "015-05" |  (substr(titre,11,16))== "015-06"
# x9= (substr(titre,11,16))== "015-07" |  (substr(titre,11,16))== "015-08" |  (substr(titre,11,16))== "015-09" |  (substr(titre,11,16))== "015-10" |  (substr(titre,11,16))== "015-11" |  (substr(titre,11,16))== "015-12"
# x17=(substr(titre,11,16))== "015-13" |  (substr(titre,11,16))== "015-14" |  (substr(titre,11,16))== "015-15" |  (substr(titre,11,16))== "015-16" |  (substr(titre,11,16))== "015-17" |  (substr(titre,11,16))== "015-18"
# x2= (substr(titre,11,16))== "169-01" |  (substr(titre,11,16))== "169-02" |  (substr(titre,11,16))== "169-03" |  (substr(titre,11,16))== "169-04" |  (substr(titre,11,16))== "169-05" |  (substr(titre,11,16))== "169-06"
# x10=(substr(titre,11,16))== "169-07" |  (substr(titre,11,16))== "169-08" |  (substr(titre,11,16))== "169-09" |  (substr(titre,11,16))== "169-10" |  (substr(titre,11,16))== "169-11" |  (substr(titre,11,16))== "169-12"
# x18=(substr(titre,11,16))== "169-13" |  (substr(titre,11,16))== "169-14" |  (substr(titre,11,16))== "169-15" |  (substr(titre,11,16))== "169-16" |  (substr(titre,11,16))== "169-17" |  (substr(titre,11,16))== "169-18"
# x3= (substr(titre,11,16))== "222-01" |  (substr(titre,11,16))== "222-02" |  (substr(titre,11,16))== "222-03" |  (substr(titre,11,16))== "222-04" |  (substr(titre,11,16))== "222-05" |  (substr(titre,11,16))== "222-06"
# x11=(substr(titre,11,16))== "222-07" |  (substr(titre,11,16))== "222-08" |  (substr(titre,11,16))== "222-09" |  (substr(titre,11,16))== "222-10" |  (substr(titre,11,16))== "222-11" |  (substr(titre,11,16))== "222-12"
# x19=(substr(titre,11,16))== "222-13" |  (substr(titre,11,16))== "222-14" |  (substr(titre,11,16))== "222-15" |  (substr(titre,11,16))== "222-16" |  (substr(titre,11,16))== "222-17" |  (substr(titre,11,16))== "222-18"
# x4= (substr(titre,11,16))== "471-01" |  (substr(titre,11,16))== "471-02" |  (substr(titre,11,16))== "471-03" |  (substr(titre,11,16))== "471-04" |  (substr(titre,11,16))== "471-05" |  (substr(titre,11,16))== "471-06"
# x12=(substr(titre,11,16))== "471-07" |  (substr(titre,11,16))== "471-08" |  (substr(titre,11,16))== "471-09" |  (substr(titre,11,16))== "471-10" |  (substr(titre,11,16))== "471-11" |  (substr(titre,11,16))== "471-12"
# x20=(substr(titre,11,16))== "471-13" |  (substr(titre,11,16))== "471-14" |  (substr(titre,11,16))== "471-15" |  (substr(titre,11,16))== "471-16" |  (substr(titre,11,16))== "471-17" |  (substr(titre,11,16))== "471-18"
# x5= (substr(titre,11,16))== "509-01" |  (substr(titre,11,16))== "509-02" |  (substr(titre,11,16))== "509-03" |  (substr(titre,11,16))== "509-04" |  (substr(titre,11,16))== "509-05" |  (substr(titre,11,16))== "509-06"
# x13=(substr(titre,11,16))== "509-07" |  (substr(titre,11,16))== "509-08" |  (substr(titre,11,16))== "509-09" |  (substr(titre,11,16))== "509-10" |  (substr(titre,11,16))== "509-11" |  (substr(titre,11,16))== "509-12"
# x21=(substr(titre,11,16))== "509-13" |  (substr(titre,11,16))== "509-14" |  (substr(titre,11,16))== "509-15" |  (substr(titre,11,16))== "509-16" |  (substr(titre,11,16))== "509-17" |  (substr(titre,11,16))== "509-18"
# x6= (substr(titre,11,16))== "525-01" |  (substr(titre,11,16))== "525-02" |  (substr(titre,11,16))== "525-03" |  (substr(titre,11,16))== "525-04" |  (substr(titre,11,16))== "525-05" |  (substr(titre,11,16))== "525-06"
# x14=(substr(titre,11,16))== "525-07" |  (substr(titre,11,16))== "525-08" |  (substr(titre,11,16))== "525-09" |  (substr(titre,11,16))== "525-10" |  (substr(titre,11,16))== "525-11" |  (substr(titre,11,16))== "525-12"
# x22=(substr(titre,11,16))== "525-13" |  (substr(titre,11,16))== "525-14" |  (substr(titre,11,16))== "525-15" |  (substr(titre,11,16))== "525-16" |  (substr(titre,11,16))== "525-17" |  (substr(titre,11,16))== "525-18"
# x7= (substr(titre,11,16))== "685-01" |  (substr(titre,11,16))== "685-02" |  (substr(titre,11,16))== "685-03" |  (substr(titre,11,16))== "685-04" |  (substr(titre,11,16))== "685-05" |  (substr(titre,11,16))== "685-06"
# x15=(substr(titre,11,16))== "685-07" |  (substr(titre,11,16))== "685-08" |  (substr(titre,11,16))== "685-09" |  (substr(titre,11,16))== "685-10" |  (substr(titre,11,16))== "685-11" |  (substr(titre,11,16))== "685-12"
# x23=(substr(titre,11,16))== "685-13" |  (substr(titre,11,16))== "685-14" |  (substr(titre,11,16))== "685-15" |  (substr(titre,11,16))== "685-16" |  (substr(titre,11,16))== "685-17" |  (substr(titre,11,16))== "685-18"
# x8= (substr(titre,11,16))== "787-01" |  (substr(titre,11,16))== "787-02" |  (substr(titre,11,16))== "787-03" |  (substr(titre,11,16))== "787-04" |  (substr(titre,11,16))== "787-05" |  (substr(titre,11,16))== "787-06"
# x16=(substr(titre,11,16))== "787-07" |  (substr(titre,11,16))== "787-08" |  (substr(titre,11,16))== "787-09" |  (substr(titre,11,16))== "787-10" |  (substr(titre,11,16))== "787-11" |  (substr(titre,11,16))== "787-12"
# x24=(substr(titre,11,16))== "787-13" |  (substr(titre,11,16))== "787-14" |  (substr(titre,11,16))== "787-15" |  (substr(titre,11,16))== "787-16" |  (substr(titre,11,16))== "787-17" |  (substr(titre,11,16))== "787-18"

# #Comparaison 471 et 877
# x1= (substr(titre,11,16))== "015-01" |  (substr(titre,11,16))== "015-02" |  (substr(titre,11,16))== "015-03" |  (substr(titre,11,16))== "015-04" |  (substr(titre,11,16))== "015-05" |  (substr(titre,11,16))== "015-06"
# x9= (substr(titre,11,16))== "015-07" |  (substr(titre,11,16))== "015-08" |  (substr(titre,11,16))== "015-09" |  (substr(titre,11,16))== "015-10" |  (substr(titre,11,16))== "015-11" |  (substr(titre,11,16))== "015-12"
# x17=(substr(titre,11,16))== "015-13" |  (substr(titre,11,16))== "015-14" |  (substr(titre,11,16))== "015-15" |  (substr(titre,11,16))== "015-16" |  (substr(titre,11,16))== "015-17" |  (substr(titre,11,16))== "015-18"
# x2= (substr(titre,11,16))== "169-01" |  (substr(titre,11,16))== "169-02" |  (substr(titre,11,16))== "169-03" |  (substr(titre,11,16))== "169-04" |  (substr(titre,11,16))== "169-05" |  (substr(titre,11,16))== "169-06"
# x10=(substr(titre,11,16))== "169-07" |  (substr(titre,11,16))== "169-08" |  (substr(titre,11,16))== "169-09" |  (substr(titre,11,16))== "169-10" |  (substr(titre,11,16))== "169-11" |  (substr(titre,11,16))== "169-12"
# x18=(substr(titre,11,16))== "169-13" |  (substr(titre,11,16))== "169-14" |  (substr(titre,11,16))== "169-15" |  (substr(titre,11,16))== "169-16" |  (substr(titre,11,16))== "169-17" |  (substr(titre,11,16))== "169-18"
# x3= (substr(titre,11,16))== "222-01" |  (substr(titre,11,16))== "222-02" |  (substr(titre,11,16))== "222-03" |  (substr(titre,11,16))== "222-04" |  (substr(titre,11,16))== "222-05" |  (substr(titre,11,16))== "222-06"
# x11=(substr(titre,11,16))== "222-07" |  (substr(titre,11,16))== "222-08" |  (substr(titre,11,16))== "222-09" |  (substr(titre,11,16))== "222-10" |  (substr(titre,11,16))== "222-11" |  (substr(titre,11,16))== "222-12"
# x19=(substr(titre,11,16))== "222-13" |  (substr(titre,11,16))== "222-14" |  (substr(titre,11,16))== "222-15" |  (substr(titre,11,16))== "222-16" |  (substr(titre,11,16))== "222-17" |  (substr(titre,11,16))== "222-18"
# x4= (substr(titre,11,16))== "471-01" |  (substr(titre,11,16))== "471-02" |  (substr(titre,11,16))== "471-03" |  (substr(titre,11,16))== "471-04" |  (substr(titre,11,16))== "471-05" |  (substr(titre,11,16))== "471-06"
# x12=(substr(titre,11,16))== "471-07" |  (substr(titre,11,16))== "471-08" |  (substr(titre,11,16))== "471-09" |  (substr(titre,11,16))== "471-10" |  (substr(titre,11,16))== "471-11" |  (substr(titre,11,16))== "471-12"
# x20=(substr(titre,11,16))== "471-13" |  (substr(titre,11,16))== "471-14" |  (substr(titre,11,16))== "471-15" |  (substr(titre,11,16))== "471-16" |  (substr(titre,11,16))== "471-17" |  (substr(titre,11,16))== "471-18"
# x5= (substr(titre,11,16))== "509-01" |  (substr(titre,11,16))== "509-02" |  (substr(titre,11,16))== "509-03" |  (substr(titre,11,16))== "509-04" |  (substr(titre,11,16))== "509-05" |  (substr(titre,11,16))== "509-06"
# x13=(substr(titre,11,16))== "509-07" |  (substr(titre,11,16))== "509-08" |  (substr(titre,11,16))== "509-09" |  (substr(titre,11,16))== "509-10" |  (substr(titre,11,16))== "509-11" |  (substr(titre,11,16))== "509-12"
# x21=(substr(titre,11,16))== "509-13" |  (substr(titre,11,16))== "509-14" |  (substr(titre,11,16))== "509-15" |  (substr(titre,11,16))== "509-16" |  (substr(titre,11,16))== "509-17" |  (substr(titre,11,16))== "509-18"
# x6= (substr(titre,11,16))== "685-01" |  (substr(titre,11,16))== "685-02" |  (substr(titre,11,16))== "685-03" |  (substr(titre,11,16))== "685-04" |  (substr(titre,11,16))== "685-05" |  (substr(titre,11,16))== "685-06"
# x14=(substr(titre,11,16))== "685-07" |  (substr(titre,11,16))== "685-08" |  (substr(titre,11,16))== "685-09" |  (substr(titre,11,16))== "685-10" |  (substr(titre,11,16))== "685-11" |  (substr(titre,11,16))== "685-12"
# x22=(substr(titre,11,16))== "685-13" |  (substr(titre,11,16))== "685-14" |  (substr(titre,11,16))== "685-15" |  (substr(titre,11,16))== "685-16" |  (substr(titre,11,16))== "685-17" |  (substr(titre,11,16))== "685-18"
# x7= (substr(titre,11,16))== "787-01" |  (substr(titre,11,16))== "787-02" |  (substr(titre,11,16))== "787-03" |  (substr(titre,11,16))== "787-04" |  (substr(titre,11,16))== "787-05" |  (substr(titre,11,16))== "787-06"
# x15=(substr(titre,11,16))== "787-07" |  (substr(titre,11,16))== "787-08" |  (substr(titre,11,16))== "787-09" |  (substr(titre,11,16))== "787-10" |  (substr(titre,11,16))== "787-11" |  (substr(titre,11,16))== "787-12"
# x23=(substr(titre,11,16))== "787-13" |  (substr(titre,11,16))== "787-14" |  (substr(titre,11,16))== "787-15" |  (substr(titre,11,16))== "787-16" |  (substr(titre,11,16))== "787-17" |  (substr(titre,11,16))== "787-18"
# x8= (substr(titre,11,16))== "877-01" |  (substr(titre,11,16))== "877-02" |  (substr(titre,11,16))== "877-03" |  (substr(titre,11,16))== "877-04" |  (substr(titre,11,16))== "877-05" |  (substr(titre,11,16))== "877-06"
# x16=(substr(titre,11,16))== "877-07" |  (substr(titre,11,16))== "877-08" |  (substr(titre,11,16))== "877-09" |  (substr(titre,11,16))== "877-10" |  (substr(titre,11,16))== "877-11" |  (substr(titre,11,16))== "877-12"
# x24=(substr(titre,11,16))== "877-13" |  (substr(titre,11,16))== "877-14" |  (substr(titre,11,16))== "877-15" |  (substr(titre,11,16))== "877-16" |  (substr(titre,11,16))== "877-17" |  (substr(titre,11,16))== "877-18"

# #Comparaison 525 et 877
# x1= (substr(titre,11,16))== "015-01" |  (substr(titre,11,16))== "015-02" |  (substr(titre,11,16))== "015-03" |  (substr(titre,11,16))== "015-04" |  (substr(titre,11,16))== "015-05" |  (substr(titre,11,16))== "015-06"
# x9= (substr(titre,11,16))== "015-07" |  (substr(titre,11,16))== "015-08" |  (substr(titre,11,16))== "015-09" |  (substr(titre,11,16))== "015-10" |  (substr(titre,11,16))== "015-11" |  (substr(titre,11,16))== "015-12"
# x17=(substr(titre,11,16))== "015-13" |  (substr(titre,11,16))== "015-14" |  (substr(titre,11,16))== "015-15" |  (substr(titre,11,16))== "015-16" |  (substr(titre,11,16))== "015-17" |  (substr(titre,11,16))== "015-18"
# x2= (substr(titre,11,16))== "169-01" |  (substr(titre,11,16))== "169-02" |  (substr(titre,11,16))== "169-03" |  (substr(titre,11,16))== "169-04" |  (substr(titre,11,16))== "169-05" |  (substr(titre,11,16))== "169-06"
# x10=(substr(titre,11,16))== "169-07" |  (substr(titre,11,16))== "169-08" |  (substr(titre,11,16))== "169-09" |  (substr(titre,11,16))== "169-10" |  (substr(titre,11,16))== "169-11" |  (substr(titre,11,16))== "169-12"
# x18=(substr(titre,11,16))== "169-13" |  (substr(titre,11,16))== "169-14" |  (substr(titre,11,16))== "169-15" |  (substr(titre,11,16))== "169-16" |  (substr(titre,11,16))== "169-17" |  (substr(titre,11,16))== "169-18"
# x3= (substr(titre,11,16))== "222-01" |  (substr(titre,11,16))== "222-02" |  (substr(titre,11,16))== "222-03" |  (substr(titre,11,16))== "222-04" |  (substr(titre,11,16))== "222-05" |  (substr(titre,11,16))== "222-06"
# x11=(substr(titre,11,16))== "222-07" |  (substr(titre,11,16))== "222-08" |  (substr(titre,11,16))== "222-09" |  (substr(titre,11,16))== "222-10" |  (substr(titre,11,16))== "222-11" |  (substr(titre,11,16))== "222-12"
# x19=(substr(titre,11,16))== "222-13" |  (substr(titre,11,16))== "222-14" |  (substr(titre,11,16))== "222-15" |  (substr(titre,11,16))== "222-16" |  (substr(titre,11,16))== "222-17" |  (substr(titre,11,16))== "222-18"
# x4= (substr(titre,11,16))== "509-01" |  (substr(titre,11,16))== "509-02" |  (substr(titre,11,16))== "509-03" |  (substr(titre,11,16))== "509-04" |  (substr(titre,11,16))== "509-05" |  (substr(titre,11,16))== "509-06"
# x12=(substr(titre,11,16))== "509-07" |  (substr(titre,11,16))== "509-08" |  (substr(titre,11,16))== "509-09" |  (substr(titre,11,16))== "509-10" |  (substr(titre,11,16))== "509-11" |  (substr(titre,11,16))== "509-12"
# x20=(substr(titre,11,16))== "509-13" |  (substr(titre,11,16))== "509-14" |  (substr(titre,11,16))== "509-15" |  (substr(titre,11,16))== "509-16" |  (substr(titre,11,16))== "509-17" |  (substr(titre,11,16))== "509-18"
# x5= (substr(titre,11,16))== "525-01" |  (substr(titre,11,16))== "525-02" |  (substr(titre,11,16))== "525-03" |  (substr(titre,11,16))== "525-04" |  (substr(titre,11,16))== "525-05" |  (substr(titre,11,16))== "525-06"
# x13=(substr(titre,11,16))== "525-07" |  (substr(titre,11,16))== "525-08" |  (substr(titre,11,16))== "525-09" |  (substr(titre,11,16))== "525-10" |  (substr(titre,11,16))== "525-11" |  (substr(titre,11,16))== "525-12"
# x21=(substr(titre,11,16))== "525-13" |  (substr(titre,11,16))== "525-14" |  (substr(titre,11,16))== "525-15" |  (substr(titre,11,16))== "525-16" |  (substr(titre,11,16))== "525-17" |  (substr(titre,11,16))== "525-18"
# x6= (substr(titre,11,16))== "685-01" |  (substr(titre,11,16))== "685-02" |  (substr(titre,11,16))== "685-03" |  (substr(titre,11,16))== "685-04" |  (substr(titre,11,16))== "685-05" |  (substr(titre,11,16))== "685-06"
# x14=(substr(titre,11,16))== "685-07" |  (substr(titre,11,16))== "685-08" |  (substr(titre,11,16))== "685-09" |  (substr(titre,11,16))== "685-10" |  (substr(titre,11,16))== "685-11" |  (substr(titre,11,16))== "685-12"
# x22=(substr(titre,11,16))== "685-13" |  (substr(titre,11,16))== "685-14" |  (substr(titre,11,16))== "685-15" |  (substr(titre,11,16))== "685-16" |  (substr(titre,11,16))== "685-17" |  (substr(titre,11,16))== "685-18"
# x7= (substr(titre,11,16))== "787-01" |  (substr(titre,11,16))== "787-02" |  (substr(titre,11,16))== "787-03" |  (substr(titre,11,16))== "787-04" |  (substr(titre,11,16))== "787-05" |  (substr(titre,11,16))== "787-06"
# x15=(substr(titre,11,16))== "787-07" |  (substr(titre,11,16))== "787-08" |  (substr(titre,11,16))== "787-09" |  (substr(titre,11,16))== "787-10" |  (substr(titre,11,16))== "787-11" |  (substr(titre,11,16))== "787-12"
# x23=(substr(titre,11,16))== "787-13" |  (substr(titre,11,16))== "787-14" |  (substr(titre,11,16))== "787-15" |  (substr(titre,11,16))== "787-16" |  (substr(titre,11,16))== "787-17" |  (substr(titre,11,16))== "787-18"
# x8= (substr(titre,11,16))== "877-01" |  (substr(titre,11,16))== "877-02" |  (substr(titre,11,16))== "877-03" |  (substr(titre,11,16))== "877-04" |  (substr(titre,11,16))== "877-05" |  (substr(titre,11,16))== "877-06"
# x16=(substr(titre,11,16))== "877-07" |  (substr(titre,11,16))== "877-08" |  (substr(titre,11,16))== "877-09" |  (substr(titre,11,16))== "877-10" |  (substr(titre,11,16))== "877-11" |  (substr(titre,11,16))== "877-12"
# x24=(substr(titre,11,16))== "877-13" |  (substr(titre,11,16))== "877-14" |  (substr(titre,11,16))== "877-15" |  (substr(titre,11,16))== "877-16" |  (substr(titre,11,16))== "877-17" |  (substr(titre,11,16))== "877-18"

# ##Comparaison 525 et 747
# x1= (substr(titre,11,16))== "015-01" |  (substr(titre,11,16))== "015-02" |  (substr(titre,11,16))== "015-03" |  (substr(titre,11,16))== "015-04" |  (substr(titre,11,16))== "015-05" |  (substr(titre,11,16))== "015-06"
# x9= (substr(titre,11,16))== "015-07" |  (substr(titre,11,16))== "015-08" |  (substr(titre,11,16))== "015-09" |  (substr(titre,11,16))== "015-10" |  (substr(titre,11,16))== "015-11" |  (substr(titre,11,16))== "015-12"
# x17=(substr(titre,11,16))== "015-13" |  (substr(titre,11,16))== "015-14" |  (substr(titre,11,16))== "015-15" |  (substr(titre,11,16))== "015-16" |  (substr(titre,11,16))== "015-17" |  (substr(titre,11,16))== "015-18"
# x2= (substr(titre,11,16))== "169-01" |  (substr(titre,11,16))== "169-02" |  (substr(titre,11,16))== "169-03" |  (substr(titre,11,16))== "169-04" |  (substr(titre,11,16))== "169-05" |  (substr(titre,11,16))== "169-06"
# x10=(substr(titre,11,16))== "169-07" |  (substr(titre,11,16))== "169-08" |  (substr(titre,11,16))== "169-09" |  (substr(titre,11,16))== "169-10" |  (substr(titre,11,16))== "169-11" |  (substr(titre,11,16))== "169-12"
# x18=(substr(titre,11,16))== "169-13" |  (substr(titre,11,16))== "169-14" |  (substr(titre,11,16))== "169-15" |  (substr(titre,11,16))== "169-16" |  (substr(titre,11,16))== "169-17" |  (substr(titre,11,16))== "169-18"
# x3= (substr(titre,11,16))== "222-01" |  (substr(titre,11,16))== "222-02" |  (substr(titre,11,16))== "222-03" |  (substr(titre,11,16))== "222-04" |  (substr(titre,11,16))== "222-05" |  (substr(titre,11,16))== "222-06"
# x11=(substr(titre,11,16))== "222-07" |  (substr(titre,11,16))== "222-08" |  (substr(titre,11,16))== "222-09" |  (substr(titre,11,16))== "222-10" |  (substr(titre,11,16))== "222-11" |  (substr(titre,11,16))== "222-12"
# x19=(substr(titre,11,16))== "222-13" |  (substr(titre,11,16))== "222-14" |  (substr(titre,11,16))== "222-15" |  (substr(titre,11,16))== "222-16" |  (substr(titre,11,16))== "222-17" |  (substr(titre,11,16))== "222-18"
# x4= (substr(titre,11,16))== "509-01" |  (substr(titre,11,16))== "509-02" |  (substr(titre,11,16))== "509-03" |  (substr(titre,11,16))== "509-04" |  (substr(titre,11,16))== "509-05" |  (substr(titre,11,16))== "509-06"
# x12=(substr(titre,11,16))== "509-07" |  (substr(titre,11,16))== "509-08" |  (substr(titre,11,16))== "509-09" |  (substr(titre,11,16))== "509-10" |  (substr(titre,11,16))== "509-11" |  (substr(titre,11,16))== "509-12"
# x20=(substr(titre,11,16))== "509-13" |  (substr(titre,11,16))== "509-14" |  (substr(titre,11,16))== "509-15" |  (substr(titre,11,16))== "509-16" |  (substr(titre,11,16))== "509-17" |  (substr(titre,11,16))== "509-18"
# x5= (substr(titre,11,16))== "525-01" |  (substr(titre,11,16))== "525-02" |  (substr(titre,11,16))== "525-03" |  (substr(titre,11,16))== "525-04" |  (substr(titre,11,16))== "525-05" |  (substr(titre,11,16))== "525-06"
# x13=(substr(titre,11,16))== "525-07" |  (substr(titre,11,16))== "525-08" |  (substr(titre,11,16))== "525-09" |  (substr(titre,11,16))== "525-10" |  (substr(titre,11,16))== "525-11" |  (substr(titre,11,16))== "525-12"
# x21=(substr(titre,11,16))== "525-13" |  (substr(titre,11,16))== "525-14" |  (substr(titre,11,16))== "525-15" |  (substr(titre,11,16))== "525-16" |  (substr(titre,11,16))== "525-17" |  (substr(titre,11,16))== "525-18"
# x6= (substr(titre,11,16))== "685-01" |  (substr(titre,11,16))== "685-02" |  (substr(titre,11,16))== "685-03" |  (substr(titre,11,16))== "685-04" |  (substr(titre,11,16))== "685-05" |  (substr(titre,11,16))== "685-06"
# x14=(substr(titre,11,16))== "685-07" |  (substr(titre,11,16))== "685-08" |  (substr(titre,11,16))== "685-09" |  (substr(titre,11,16))== "685-10" |  (substr(titre,11,16))== "685-11" |  (substr(titre,11,16))== "685-12"
# x22=(substr(titre,11,16))== "685-13" |  (substr(titre,11,16))== "685-14" |  (substr(titre,11,16))== "685-15" |  (substr(titre,11,16))== "685-16" |  (substr(titre,11,16))== "685-17" |  (substr(titre,11,16))== "685-18"
# x7= (substr(titre,11,16))== "747-01" |  (substr(titre,11,16))== "747-02" |  (substr(titre,11,16))== "747-03" |  (substr(titre,11,16))== "747-04" |  (substr(titre,11,16))== "747-05" |  (substr(titre,11,16))== "747-06"
# x15=(substr(titre,11,16))== "747-07" |  (substr(titre,11,16))== "747-08" |  (substr(titre,11,16))== "747-09" |  (substr(titre,11,16))== "747-10" |  (substr(titre,11,16))== "747-11" |  (substr(titre,11,16))== "747-12"
# x23=(substr(titre,11,16))== "747-13" |  (substr(titre,11,16))== "747-14" |  (substr(titre,11,16))== "747-15" |  (substr(titre,11,16))== "747-16" |  (substr(titre,11,16))== "747-17" |  (substr(titre,11,16))== "747-18"
# x8= (substr(titre,11,16))== "787-01" |  (substr(titre,11,16))== "787-02" |  (substr(titre,11,16))== "787-03" |  (substr(titre,11,16))== "787-04" |  (substr(titre,11,16))== "787-05" |  (substr(titre,11,16))== "787-06"
# x16=(substr(titre,11,16))== "787-07" |  (substr(titre,11,16))== "787-08" |  (substr(titre,11,16))== "787-09" |  (substr(titre,11,16))== "787-10" |  (substr(titre,11,16))== "787-11" |  (substr(titre,11,16))== "787-12"
# x24=(substr(titre,11,16))== "787-13" |  (substr(titre,11,16))== "787-14" |  (substr(titre,11,16))== "787-15" |  (substr(titre,11,16))== "787-16" |  (substr(titre,11,16))== "787-17" |  (substr(titre,11,16))== "787-18"

# #Comparaison 471 et 747
# x1= (substr(titre,11,16))== "015-01" |  (substr(titre,11,16))== "015-02" |  (substr(titre,11,16))== "015-03" |  (substr(titre,11,16))== "015-04" |  (substr(titre,11,16))== "015-05" |  (substr(titre,11,16))== "015-06"
# x9= (substr(titre,11,16))== "015-07" |  (substr(titre,11,16))== "015-08" |  (substr(titre,11,16))== "015-09" |  (substr(titre,11,16))== "015-10" |  (substr(titre,11,16))== "015-11" |  (substr(titre,11,16))== "015-12"
# x17=(substr(titre,11,16))== "015-13" |  (substr(titre,11,16))== "015-14" |  (substr(titre,11,16))== "015-15" |  (substr(titre,11,16))== "015-16" |  (substr(titre,11,16))== "015-17" |  (substr(titre,11,16))== "015-18"
# x2= (substr(titre,11,16))== "169-01" |  (substr(titre,11,16))== "169-02" |  (substr(titre,11,16))== "169-03" |  (substr(titre,11,16))== "169-04" |  (substr(titre,11,16))== "169-05" |  (substr(titre,11,16))== "169-06"
# x10=(substr(titre,11,16))== "169-07" |  (substr(titre,11,16))== "169-08" |  (substr(titre,11,16))== "169-09" |  (substr(titre,11,16))== "169-10" |  (substr(titre,11,16))== "169-11" |  (substr(titre,11,16))== "169-12"
# x18=(substr(titre,11,16))== "169-13" |  (substr(titre,11,16))== "169-14" |  (substr(titre,11,16))== "169-15" |  (substr(titre,11,16))== "169-16" |  (substr(titre,11,16))== "169-17" |  (substr(titre,11,16))== "169-18"
# x3= (substr(titre,11,16))== "222-01" |  (substr(titre,11,16))== "222-02" |  (substr(titre,11,16))== "222-03" |  (substr(titre,11,16))== "222-04" |  (substr(titre,11,16))== "222-05" |  (substr(titre,11,16))== "222-06"
# x11=(substr(titre,11,16))== "222-07" |  (substr(titre,11,16))== "222-08" |  (substr(titre,11,16))== "222-09" |  (substr(titre,11,16))== "222-10" |  (substr(titre,11,16))== "222-11" |  (substr(titre,11,16))== "222-12"
# x19=(substr(titre,11,16))== "222-13" |  (substr(titre,11,16))== "222-14" |  (substr(titre,11,16))== "222-15" |  (substr(titre,11,16))== "222-16" |  (substr(titre,11,16))== "222-17" |  (substr(titre,11,16))== "222-18"
# x4= (substr(titre,11,16))== "471-01" |  (substr(titre,11,16))== "471-02" |  (substr(titre,11,16))== "471-03" |  (substr(titre,11,16))== "471-04" |  (substr(titre,11,16))== "471-05" |  (substr(titre,11,16))== "471-06"
# x12=(substr(titre,11,16))== "471-07" |  (substr(titre,11,16))== "471-08" |  (substr(titre,11,16))== "471-09" |  (substr(titre,11,16))== "471-10" |  (substr(titre,11,16))== "471-11" |  (substr(titre,11,16))== "471-12"
# x20=(substr(titre,11,16))== "471-13" |  (substr(titre,11,16))== "471-14" |  (substr(titre,11,16))== "471-15" |  (substr(titre,11,16))== "471-16" |  (substr(titre,11,16))== "471-17" |  (substr(titre,11,16))== "471-18"
# x5= (substr(titre,11,16))== "509-01" |  (substr(titre,11,16))== "509-02" |  (substr(titre,11,16))== "509-03" |  (substr(titre,11,16))== "509-04" |  (substr(titre,11,16))== "509-05" |  (substr(titre,11,16))== "509-06"
# x13=(substr(titre,11,16))== "509-07" |  (substr(titre,11,16))== "509-08" |  (substr(titre,11,16))== "509-09" |  (substr(titre,11,16))== "509-10" |  (substr(titre,11,16))== "509-11" |  (substr(titre,11,16))== "509-12"
# x21=(substr(titre,11,16))== "509-13" |  (substr(titre,11,16))== "509-14" |  (substr(titre,11,16))== "509-15" |  (substr(titre,11,16))== "509-16" |  (substr(titre,11,16))== "509-17" |  (substr(titre,11,16))== "509-18"
# x6= (substr(titre,11,16))== "685-01" |  (substr(titre,11,16))== "685-02" |  (substr(titre,11,16))== "685-03" |  (substr(titre,11,16))== "685-04" |  (substr(titre,11,16))== "685-05" |  (substr(titre,11,16))== "685-06"
# x14=(substr(titre,11,16))== "685-07" |  (substr(titre,11,16))== "685-08" |  (substr(titre,11,16))== "685-09" |  (substr(titre,11,16))== "685-10" |  (substr(titre,11,16))== "685-11" |  (substr(titre,11,16))== "685-12"
# x22=(substr(titre,11,16))== "685-13" |  (substr(titre,11,16))== "685-14" |  (substr(titre,11,16))== "685-15" |  (substr(titre,11,16))== "685-16" |  (substr(titre,11,16))== "685-17" |  (substr(titre,11,16))== "685-18"
# x7= (substr(titre,11,16))== "747-01" |  (substr(titre,11,16))== "747-02" |  (substr(titre,11,16))== "747-03" |  (substr(titre,11,16))== "747-04" |  (substr(titre,11,16))== "747-05" |  (substr(titre,11,16))== "747-06"
# x15=(substr(titre,11,16))== "747-07" |  (substr(titre,11,16))== "747-08" |  (substr(titre,11,16))== "747-09" |  (substr(titre,11,16))== "747-10" |  (substr(titre,11,16))== "747-11" |  (substr(titre,11,16))== "747-12"
# x23=(substr(titre,11,16))== "747-13" |  (substr(titre,11,16))== "747-14" |  (substr(titre,11,16))== "747-15" |  (substr(titre,11,16))== "747-16" |  (substr(titre,11,16))== "747-17" |  (substr(titre,11,16))== "747-18"
# x8= (substr(titre,11,16))== "787-01" |  (substr(titre,11,16))== "787-02" |  (substr(titre,11,16))== "787-03" |  (substr(titre,11,16))== "787-04" |  (substr(titre,11,16))== "787-05" |  (substr(titre,11,16))== "787-06"
# x16=(substr(titre,11,16))== "787-07" |  (substr(titre,11,16))== "787-08" |  (substr(titre,11,16))== "787-09" |  (substr(titre,11,16))== "787-10" |  (substr(titre,11,16))== "787-11" |  (substr(titre,11,16))== "787-12"
# x24=(substr(titre,11,16))== "787-13" |  (substr(titre,11,16))== "787-14" |  (substr(titre,11,16))== "787-15" |  (substr(titre,11,16))== "787-16" |  (substr(titre,11,16))== "787-17" |  (substr(titre,11,16))== "787-18"

# #Comparaison 015 et 169
# x1= (substr(titre,11,16))== "015-01" |  (substr(titre,11,16))== "015-02" |  (substr(titre,11,16))== "015-03" |  (substr(titre,11,16))== "015-04" |  (substr(titre,11,16))== "015-05" |  (substr(titre,11,16))== "015-06"
# x10=(substr(titre,11,16))== "015-07" |  (substr(titre,11,16))== "015-08" |  (substr(titre,11,16))== "015-09" |  (substr(titre,11,16))== "015-10" |  (substr(titre,11,16))== "015-11" |  (substr(titre,11,16))== "015-12"
# x19=(substr(titre,11,16))== "015-13" |  (substr(titre,11,16))== "015-14" |  (substr(titre,11,16))== "015-15" |  (substr(titre,11,16))== "015-16" |  (substr(titre,11,16))== "015-17" |  (substr(titre,11,16))== "015-18"
# x2= (substr(titre,11,16))== "169-01" |  (substr(titre,11,16))== "169-02" |  (substr(titre,11,16))== "169-03" |  (substr(titre,11,16))== "169-04" |  (substr(titre,11,16))== "169-05" |  (substr(titre,11,16))== "169-06"
# x11=(substr(titre,11,16))== "169-07" |  (substr(titre,11,16))== "169-08" |  (substr(titre,11,16))== "169-09" |  (substr(titre,11,16))== "169-10" |  (substr(titre,11,16))== "169-11" |  (substr(titre,11,16))== "169-12"
# x20=(substr(titre,11,16))== "169-13" |  (substr(titre,11,16))== "169-14" |  (substr(titre,11,16))== "169-15" |  (substr(titre,11,16))== "169-16" |  (substr(titre,11,16))== "169-17" |  (substr(titre,11,16))== "169-18"
# x3= (substr(titre,11,16))== "222-01" |  (substr(titre,11,16))== "222-02" |  (substr(titre,11,16))== "222-03" |  (substr(titre,11,16))== "222-04" |  (substr(titre,11,16))== "222-05" |  (substr(titre,11,16))== "222-06"
# x12=(substr(titre,11,16))== "222-07" |  (substr(titre,11,16))== "222-08" |  (substr(titre,11,16))== "222-09" |  (substr(titre,11,16))== "222-10" |  (substr(titre,11,16))== "222-11" |  (substr(titre,11,16))== "222-12"
# x21=(substr(titre,11,16))== "222-13" |  (substr(titre,11,16))== "222-14" |  (substr(titre,11,16))== "222-15" |  (substr(titre,11,16))== "222-16" |  (substr(titre,11,16))== "222-17" |  (substr(titre,11,16))== "222-18"
# x4= (substr(titre,11,16))== "471-01" |  (substr(titre,11,16))== "471-02" |  (substr(titre,11,16))== "471-03" |  (substr(titre,11,16))== "471-04" |  (substr(titre,11,16))== "471-05" |  (substr(titre,11,16))== "471-06"
# x13=(substr(titre,11,16))== "471-07" |  (substr(titre,11,16))== "471-08" |  (substr(titre,11,16))== "471-09" |  (substr(titre,11,16))== "471-10" |  (substr(titre,11,16))== "471-11" |  (substr(titre,11,16))== "471-12"
# x22=(substr(titre,11,16))== "471-13" |  (substr(titre,11,16))== "471-14" |  (substr(titre,11,16))== "471-15" |  (substr(titre,11,16))== "471-16" |  (substr(titre,11,16))== "471-17" |  (substr(titre,11,16))== "471-18"
# x5= (substr(titre,11,16))== "509-01" |  (substr(titre,11,16))== "509-02" |  (substr(titre,11,16))== "509-03" |  (substr(titre,11,16))== "509-04" |  (substr(titre,11,16))== "509-05" |  (substr(titre,11,16))== "509-06"
# x14=(substr(titre,11,16))== "509-07" |  (substr(titre,11,16))== "509-08" |  (substr(titre,11,16))== "509-09" |  (substr(titre,11,16))== "509-10" |  (substr(titre,11,16))== "509-11" |  (substr(titre,11,16))== "509-12"
# x23=(substr(titre,11,16))== "509-13" |  (substr(titre,11,16))== "509-14" |  (substr(titre,11,16))== "509-15" |  (substr(titre,11,16))== "509-16" |  (substr(titre,11,16))== "509-17" |  (substr(titre,11,16))== "509-18"
# x6= (substr(titre,11,16))== "525-01" |  (substr(titre,11,16))== "525-02" |  (substr(titre,11,16))== "525-03" |  (substr(titre,11,16))== "525-04" |  (substr(titre,11,16))== "525-05" |  (substr(titre,11,16))== "525-06"
# x15=(substr(titre,11,16))== "525-07" |  (substr(titre,11,16))== "525-08" |  (substr(titre,11,16))== "525-09" |  (substr(titre,11,16))== "525-10" |  (substr(titre,11,16))== "525-11" |  (substr(titre,11,16))== "525-12"
# x24=(substr(titre,11,16))== "525-13" |  (substr(titre,11,16))== "525-14" |  (substr(titre,11,16))== "525-15" |  (substr(titre,11,16))== "525-16" |  (substr(titre,11,16))== "525-17" |  (substr(titre,11,16))== "525-18"
# x7= (substr(titre,11,16))== "747-01" |  (substr(titre,11,16))== "747-02" |  (substr(titre,11,16))== "747-03" |  (substr(titre,11,16))== "747-04" |  (substr(titre,11,16))== "747-05" |  (substr(titre,11,16))== "747-06"
# x16=(substr(titre,11,16))== "747-07" |  (substr(titre,11,16))== "747-08" |  (substr(titre,11,16))== "747-09" |  (substr(titre,11,16))== "747-10" |  (substr(titre,11,16))== "747-11" |  (substr(titre,11,16))== "747-12"
# x25=(substr(titre,11,16))== "747-13" |  (substr(titre,11,16))== "747-14" |  (substr(titre,11,16))== "747-15" |  (substr(titre,11,16))== "747-16" |  (substr(titre,11,16))== "747-17" |  (substr(titre,11,16))== "747-18"
# x8= (substr(titre,11,16))== "787-01" |  (substr(titre,11,16))== "787-02" |  (substr(titre,11,16))== "787-03" |  (substr(titre,11,16))== "787-04" |  (substr(titre,11,16))== "787-05" |  (substr(titre,11,16))== "787-06"
# x17=(substr(titre,11,16))== "787-07" |  (substr(titre,11,16))== "787-08" |  (substr(titre,11,16))== "787-09" |  (substr(titre,11,16))== "787-10" |  (substr(titre,11,16))== "787-11" |  (substr(titre,11,16))== "787-12"
# x26=(substr(titre,11,16))== "787-13" |  (substr(titre,11,16))== "787-14" |  (substr(titre,11,16))== "787-15" |  (substr(titre,11,16))== "787-16" |  (substr(titre,11,16))== "787-17" |  (substr(titre,11,16))== "787-18"
# x9= (substr(titre,11,16))== "877-01" |  (substr(titre,11,16))== "877-02" |  (substr(titre,11,16))== "877-03" |  (substr(titre,11,16))== "877-04" |  (substr(titre,11,16))== "877-05" |  (substr(titre,11,16))== "877-06"
# x18=(substr(titre,11,16))== "877-07" |  (substr(titre,11,16))== "877-08" |  (substr(titre,11,16))== "877-09" |  (substr(titre,11,16))== "877-10" |  (substr(titre,11,16))== "877-11" |  (substr(titre,11,16))== "877-12"
# x27=(substr(titre,11,16))== "877-13" |  (substr(titre,11,16))== "877-14" |  (substr(titre,11,16))== "877-15" |  (substr(titre,11,16))== "877-16" |  (substr(titre,11,16))== "877-17" |  (substr(titre,11,16))== "877-18"

# #Comparaison 015 et 685
# x1= (substr(titre,11,16))== "015-01" |  (substr(titre,11,16))== "015-02" |  (substr(titre,11,16))== "015-03" |  (substr(titre,11,16))== "015-04" |  (substr(titre,11,16))== "015-05" |  (substr(titre,11,16))== "015-06"
# x10=(substr(titre,11,16))== "015-07" |  (substr(titre,11,16))== "015-08" |  (substr(titre,11,16))== "015-09" |  (substr(titre,11,16))== "015-10" |  (substr(titre,11,16))== "015-11" |  (substr(titre,11,16))== "015-12"
# x19=(substr(titre,11,16))== "015-13" |  (substr(titre,11,16))== "015-14" |  (substr(titre,11,16))== "015-15" |  (substr(titre,11,16))== "015-16" |  (substr(titre,11,16))== "015-17" |  (substr(titre,11,16))== "015-18"
# x2= (substr(titre,11,16))== "222-01" |  (substr(titre,11,16))== "222-02" |  (substr(titre,11,16))== "222-03" |  (substr(titre,11,16))== "222-04" |  (substr(titre,11,16))== "222-05" |  (substr(titre,11,16))== "222-06"
# x11=(substr(titre,11,16))== "222-07" |  (substr(titre,11,16))== "222-08" |  (substr(titre,11,16))== "222-09" |  (substr(titre,11,16))== "222-10" |  (substr(titre,11,16))== "222-11" |  (substr(titre,11,16))== "222-12"
# x20=(substr(titre,11,16))== "222-13" |  (substr(titre,11,16))== "222-14" |  (substr(titre,11,16))== "222-15" |  (substr(titre,11,16))== "222-16" |  (substr(titre,11,16))== "222-17" |  (substr(titre,11,16))== "222-18"
# x3= (substr(titre,11,16))== "471-01" |  (substr(titre,11,16))== "471-02" |  (substr(titre,11,16))== "471-03" |  (substr(titre,11,16))== "471-04" |  (substr(titre,11,16))== "471-05" |  (substr(titre,11,16))== "471-06"
# x12=(substr(titre,11,16))== "471-07" |  (substr(titre,11,16))== "471-08" |  (substr(titre,11,16))== "471-09" |  (substr(titre,11,16))== "471-10" |  (substr(titre,11,16))== "471-11" |  (substr(titre,11,16))== "471-12"
# x21=(substr(titre,11,16))== "471-13" |  (substr(titre,11,16))== "471-14" |  (substr(titre,11,16))== "471-15" |  (substr(titre,11,16))== "471-16" |  (substr(titre,11,16))== "471-17" |  (substr(titre,11,16))== "471-18"
# x4= (substr(titre,11,16))== "509-01" |  (substr(titre,11,16))== "509-02" |  (substr(titre,11,16))== "509-03" |  (substr(titre,11,16))== "509-04" |  (substr(titre,11,16))== "509-05" |  (substr(titre,11,16))== "509-06"
# x13=(substr(titre,11,16))== "509-07" |  (substr(titre,11,16))== "509-08" |  (substr(titre,11,16))== "509-09" |  (substr(titre,11,16))== "509-10" |  (substr(titre,11,16))== "509-11" |  (substr(titre,11,16))== "509-12"
# x22=(substr(titre,11,16))== "509-13" |  (substr(titre,11,16))== "509-14" |  (substr(titre,11,16))== "509-15" |  (substr(titre,11,16))== "509-16" |  (substr(titre,11,16))== "509-17" |  (substr(titre,11,16))== "509-18"
# x5= (substr(titre,11,16))== "525-01" |  (substr(titre,11,16))== "525-02" |  (substr(titre,11,16))== "525-03" |  (substr(titre,11,16))== "525-04" |  (substr(titre,11,16))== "525-05" |  (substr(titre,11,16))== "525-06"
# x14=(substr(titre,11,16))== "525-07" |  (substr(titre,11,16))== "525-08" |  (substr(titre,11,16))== "525-09" |  (substr(titre,11,16))== "525-10" |  (substr(titre,11,16))== "525-11" |  (substr(titre,11,16))== "525-12"
# x23=(substr(titre,11,16))== "525-13" |  (substr(titre,11,16))== "525-14" |  (substr(titre,11,16))== "525-15" |  (substr(titre,11,16))== "525-16" |  (substr(titre,11,16))== "525-17" |  (substr(titre,11,16))== "525-18"
# x6= (substr(titre,11,16))== "685-01" |  (substr(titre,11,16))== "685-02" |  (substr(titre,11,16))== "685-03" |  (substr(titre,11,16))== "685-04" |  (substr(titre,11,16))== "685-05" |  (substr(titre,11,16))== "685-06"
# x15=(substr(titre,11,16))== "685-07" |  (substr(titre,11,16))== "685-08" |  (substr(titre,11,16))== "685-09" |  (substr(titre,11,16))== "685-10" |  (substr(titre,11,16))== "685-11" |  (substr(titre,11,16))== "685-12"
# x24=(substr(titre,11,16))== "685-13" |  (substr(titre,11,16))== "685-14" |  (substr(titre,11,16))== "685-15" |  (substr(titre,11,16))== "685-16" |  (substr(titre,11,16))== "685-17" |  (substr(titre,11,16))== "685-18"
# x7= (substr(titre,11,16))== "747-01" |  (substr(titre,11,16))== "747-02" |  (substr(titre,11,16))== "747-03" |  (substr(titre,11,16))== "747-04" |  (substr(titre,11,16))== "747-05" |  (substr(titre,11,16))== "747-06"
# x16=(substr(titre,11,16))== "747-07" |  (substr(titre,11,16))== "747-08" |  (substr(titre,11,16))== "747-09" |  (substr(titre,11,16))== "747-10" |  (substr(titre,11,16))== "747-11" |  (substr(titre,11,16))== "747-12"
# x25=(substr(titre,11,16))== "747-13" |  (substr(titre,11,16))== "747-14" |  (substr(titre,11,16))== "747-15" |  (substr(titre,11,16))== "747-16" |  (substr(titre,11,16))== "747-17" |  (substr(titre,11,16))== "747-18"
# x8= (substr(titre,11,16))== "787-01" |  (substr(titre,11,16))== "787-02" |  (substr(titre,11,16))== "787-03" |  (substr(titre,11,16))== "787-04" |  (substr(titre,11,16))== "787-05" |  (substr(titre,11,16))== "787-06"
# x17=(substr(titre,11,16))== "787-07" |  (substr(titre,11,16))== "787-08" |  (substr(titre,11,16))== "787-09" |  (substr(titre,11,16))== "787-10" |  (substr(titre,11,16))== "787-11" |  (substr(titre,11,16))== "787-12"
# x26=(substr(titre,11,16))== "787-13" |  (substr(titre,11,16))== "787-14" |  (substr(titre,11,16))== "787-15" |  (substr(titre,11,16))== "787-16" |  (substr(titre,11,16))== "787-17" |  (substr(titre,11,16))== "787-18"
# x9= (substr(titre,11,16))== "877-01" |  (substr(titre,11,16))== "877-02" |  (substr(titre,11,16))== "877-03" |  (substr(titre,11,16))== "877-04" |  (substr(titre,11,16))== "877-05" |  (substr(titre,11,16))== "877-06"
# x18=(substr(titre,11,16))== "877-07" |  (substr(titre,11,16))== "877-08" |  (substr(titre,11,16))== "877-09" |  (substr(titre,11,16))== "877-10" |  (substr(titre,11,16))== "877-11" |  (substr(titre,11,16))== "877-12"
# x27=(substr(titre,11,16))== "877-13" |  (substr(titre,11,16))== "877-14" |  (substr(titre,11,16))== "877-15" |  (substr(titre,11,16))== "877-16" |  (substr(titre,11,16))== "877-17" |  (substr(titre,11,16))== "877-18"

# #Comparaison 169 et 685
# x1= (substr(titre,11,16))== "169-01" |  (substr(titre,11,16))== "169-02" |  (substr(titre,11,16))== "169-03" |  (substr(titre,11,16))== "169-04" |  (substr(titre,11,16))== "169-05" |  (substr(titre,11,16))== "169-06"
# x10=(substr(titre,11,16))== "169-07" |  (substr(titre,11,16))== "169-08" |  (substr(titre,11,16))== "169-09" |  (substr(titre,11,16))== "169-10" |  (substr(titre,11,16))== "169-11" |  (substr(titre,11,16))== "169-12"
# x19=(substr(titre,11,16))== "169-13" |  (substr(titre,11,16))== "169-14" |  (substr(titre,11,16))== "169-15" |  (substr(titre,11,16))== "169-16" |  (substr(titre,11,16))== "169-17" |  (substr(titre,11,16))== "169-18"
# x2= (substr(titre,11,16))== "222-01" |  (substr(titre,11,16))== "222-02" |  (substr(titre,11,16))== "222-03" |  (substr(titre,11,16))== "222-04" |  (substr(titre,11,16))== "222-05" |  (substr(titre,11,16))== "222-06"
# x11=(substr(titre,11,16))== "222-07" |  (substr(titre,11,16))== "222-08" |  (substr(titre,11,16))== "222-09" |  (substr(titre,11,16))== "222-10" |  (substr(titre,11,16))== "222-11" |  (substr(titre,11,16))== "222-12"
# x20=(substr(titre,11,16))== "222-13" |  (substr(titre,11,16))== "222-14" |  (substr(titre,11,16))== "222-15" |  (substr(titre,11,16))== "222-16" |  (substr(titre,11,16))== "222-17" |  (substr(titre,11,16))== "222-18"
# x3= (substr(titre,11,16))== "471-01" |  (substr(titre,11,16))== "471-02" |  (substr(titre,11,16))== "471-03" |  (substr(titre,11,16))== "471-04" |  (substr(titre,11,16))== "471-05" |  (substr(titre,11,16))== "471-06"
# x12=(substr(titre,11,16))== "471-07" |  (substr(titre,11,16))== "471-08" |  (substr(titre,11,16))== "471-09" |  (substr(titre,11,16))== "471-10" |  (substr(titre,11,16))== "471-11" |  (substr(titre,11,16))== "471-12"
# x21=(substr(titre,11,16))== "471-13" |  (substr(titre,11,16))== "471-14" |  (substr(titre,11,16))== "471-15" |  (substr(titre,11,16))== "471-16" |  (substr(titre,11,16))== "471-17" |  (substr(titre,11,16))== "471-18"
# x4= (substr(titre,11,16))== "509-01" |  (substr(titre,11,16))== "509-02" |  (substr(titre,11,16))== "509-03" |  (substr(titre,11,16))== "509-04" |  (substr(titre,11,16))== "509-05" |  (substr(titre,11,16))== "509-06"
# x13=(substr(titre,11,16))== "509-07" |  (substr(titre,11,16))== "509-08" |  (substr(titre,11,16))== "509-09" |  (substr(titre,11,16))== "509-10" |  (substr(titre,11,16))== "509-11" |  (substr(titre,11,16))== "509-12"
# x22=(substr(titre,11,16))== "509-13" |  (substr(titre,11,16))== "509-14" |  (substr(titre,11,16))== "509-15" |  (substr(titre,11,16))== "509-16" |  (substr(titre,11,16))== "509-17" |  (substr(titre,11,16))== "509-18"
# x5= (substr(titre,11,16))== "525-01" |  (substr(titre,11,16))== "525-02" |  (substr(titre,11,16))== "525-03" |  (substr(titre,11,16))== "525-04" |  (substr(titre,11,16))== "525-05" |  (substr(titre,11,16))== "525-06"
# x14=(substr(titre,11,16))== "525-07" |  (substr(titre,11,16))== "525-08" |  (substr(titre,11,16))== "525-09" |  (substr(titre,11,16))== "525-10" |  (substr(titre,11,16))== "525-11" |  (substr(titre,11,16))== "525-12"
# x23=(substr(titre,11,16))== "525-13" |  (substr(titre,11,16))== "525-14" |  (substr(titre,11,16))== "525-15" |  (substr(titre,11,16))== "525-16" |  (substr(titre,11,16))== "525-17" |  (substr(titre,11,16))== "525-18"
# x6= (substr(titre,11,16))== "685-01" |  (substr(titre,11,16))== "685-02" |  (substr(titre,11,16))== "685-03" |  (substr(titre,11,16))== "685-04" |  (substr(titre,11,16))== "685-05" |  (substr(titre,11,16))== "685-06"
# x15=(substr(titre,11,16))== "685-07" |  (substr(titre,11,16))== "685-08" |  (substr(titre,11,16))== "685-09" |  (substr(titre,11,16))== "685-10" |  (substr(titre,11,16))== "685-11" |  (substr(titre,11,16))== "685-12"
# x24=(substr(titre,11,16))== "685-13" |  (substr(titre,11,16))== "685-14" |  (substr(titre,11,16))== "685-15" |  (substr(titre,11,16))== "685-16" |  (substr(titre,11,16))== "685-17" |  (substr(titre,11,16))== "685-18"
# x7= (substr(titre,11,16))== "747-01" |  (substr(titre,11,16))== "747-02" |  (substr(titre,11,16))== "747-03" |  (substr(titre,11,16))== "747-04" |  (substr(titre,11,16))== "747-05" |  (substr(titre,11,16))== "747-06"
# x16=(substr(titre,11,16))== "747-07" |  (substr(titre,11,16))== "747-08" |  (substr(titre,11,16))== "747-09" |  (substr(titre,11,16))== "747-10" |  (substr(titre,11,16))== "747-11" |  (substr(titre,11,16))== "747-12"
# x25=(substr(titre,11,16))== "747-13" |  (substr(titre,11,16))== "747-14" |  (substr(titre,11,16))== "747-15" |  (substr(titre,11,16))== "747-16" |  (substr(titre,11,16))== "747-17" |  (substr(titre,11,16))== "747-18"
# x8= (substr(titre,11,16))== "787-01" |  (substr(titre,11,16))== "787-02" |  (substr(titre,11,16))== "787-03" |  (substr(titre,11,16))== "787-04" |  (substr(titre,11,16))== "787-05" |  (substr(titre,11,16))== "787-06"
# x17=(substr(titre,11,16))== "787-07" |  (substr(titre,11,16))== "787-08" |  (substr(titre,11,16))== "787-09" |  (substr(titre,11,16))== "787-10" |  (substr(titre,11,16))== "787-11" |  (substr(titre,11,16))== "787-12"
# x26=(substr(titre,11,16))== "787-13" |  (substr(titre,11,16))== "787-14" |  (substr(titre,11,16))== "787-15" |  (substr(titre,11,16))== "787-16" |  (substr(titre,11,16))== "787-17" |  (substr(titre,11,16))== "787-18"
# x9= (substr(titre,11,16))== "877-01" |  (substr(titre,11,16))== "877-02" |  (substr(titre,11,16))== "877-03" |  (substr(titre,11,16))== "877-04" |  (substr(titre,11,16))== "877-05" |  (substr(titre,11,16))== "877-06"
# x18=(substr(titre,11,16))== "877-07" |  (substr(titre,11,16))== "877-08" |  (substr(titre,11,16))== "877-09" |  (substr(titre,11,16))== "877-10" |  (substr(titre,11,16))== "877-11" |  (substr(titre,11,16))== "877-12"
# x27=(substr(titre,11,16))== "877-13" |  (substr(titre,11,16))== "877-14" |  (substr(titre,11,16))== "877-15" |  (substr(titre,11,16))== "877-16" |  (substr(titre,11,16))== "877-17" |  (substr(titre,11,16))== "877-18"

# #Comparaison 222 et 509
# x1= (substr(titre,11,16))== "015-01" |  (substr(titre,11,16))== "015-02" |  (substr(titre,11,16))== "015-03" |  (substr(titre,11,16))== "015-04" |  (substr(titre,11,16))== "015-05" |  (substr(titre,11,16))== "015-06"
# x10=(substr(titre,11,16))== "015-07" |  (substr(titre,11,16))== "015-08" |  (substr(titre,11,16))== "015-09" |  (substr(titre,11,16))== "015-10" |  (substr(titre,11,16))== "015-11" |  (substr(titre,11,16))== "015-12"
# x19=(substr(titre,11,16))== "015-13" |  (substr(titre,11,16))== "015-14" |  (substr(titre,11,16))== "015-15" |  (substr(titre,11,16))== "015-16" |  (substr(titre,11,16))== "015-17" |  (substr(titre,11,16))== "015-18"
# x2= (substr(titre,11,16))== "169-01" |  (substr(titre,11,16))== "169-02" |  (substr(titre,11,16))== "169-03" |  (substr(titre,11,16))== "169-04" |  (substr(titre,11,16))== "169-05" |  (substr(titre,11,16))== "169-06"
# x11=(substr(titre,11,16))== "169-07" |  (substr(titre,11,16))== "169-08" |  (substr(titre,11,16))== "169-09" |  (substr(titre,11,16))== "169-10" |  (substr(titre,11,16))== "169-11" |  (substr(titre,11,16))== "169-12"
# x20=(substr(titre,11,16))== "169-13" |  (substr(titre,11,16))== "169-14" |  (substr(titre,11,16))== "169-15" |  (substr(titre,11,16))== "169-16" |  (substr(titre,11,16))== "169-17" |  (substr(titre,11,16))== "169-18"
# x3= (substr(titre,11,16))== "222-01" |  (substr(titre,11,16))== "222-02" |  (substr(titre,11,16))== "222-03" |  (substr(titre,11,16))== "222-04" |  (substr(titre,11,16))== "222-05" |  (substr(titre,11,16))== "222-06"
# x12=(substr(titre,11,16))== "222-07" |  (substr(titre,11,16))== "222-08" |  (substr(titre,11,16))== "222-09" |  (substr(titre,11,16))== "222-10" |  (substr(titre,11,16))== "222-11" |  (substr(titre,11,16))== "222-12"
# x21=(substr(titre,11,16))== "222-13" |  (substr(titre,11,16))== "222-14" |  (substr(titre,11,16))== "222-15" |  (substr(titre,11,16))== "222-16" |  (substr(titre,11,16))== "222-17" |  (substr(titre,11,16))== "222-18"
# x4= (substr(titre,11,16))== "471-01" |  (substr(titre,11,16))== "471-02" |  (substr(titre,11,16))== "471-03" |  (substr(titre,11,16))== "471-04" |  (substr(titre,11,16))== "471-05" |  (substr(titre,11,16))== "471-06"
# x13=(substr(titre,11,16))== "471-07" |  (substr(titre,11,16))== "471-08" |  (substr(titre,11,16))== "471-09" |  (substr(titre,11,16))== "471-10" |  (substr(titre,11,16))== "471-11" |  (substr(titre,11,16))== "471-12"
# x22=(substr(titre,11,16))== "471-13" |  (substr(titre,11,16))== "471-14" |  (substr(titre,11,16))== "471-15" |  (substr(titre,11,16))== "471-16" |  (substr(titre,11,16))== "471-17" |  (substr(titre,11,16))== "471-18"
# x5= (substr(titre,11,16))== "509-01" |  (substr(titre,11,16))== "509-02" |  (substr(titre,11,16))== "509-03" |  (substr(titre,11,16))== "509-04" |  (substr(titre,11,16))== "509-05" |  (substr(titre,11,16))== "509-06"
# x14=(substr(titre,11,16))== "509-07" |  (substr(titre,11,16))== "509-08" |  (substr(titre,11,16))== "509-09" |  (substr(titre,11,16))== "509-10" |  (substr(titre,11,16))== "509-11" |  (substr(titre,11,16))== "509-12"
# x23=(substr(titre,11,16))== "509-13" |  (substr(titre,11,16))== "509-14" |  (substr(titre,11,16))== "509-15" |  (substr(titre,11,16))== "509-16" |  (substr(titre,11,16))== "509-17" |  (substr(titre,11,16))== "509-18"
# x6= (substr(titre,11,16))== "525-01" |  (substr(titre,11,16))== "525-02" |  (substr(titre,11,16))== "525-03" |  (substr(titre,11,16))== "525-04" |  (substr(titre,11,16))== "525-05" |  (substr(titre,11,16))== "525-06"
# x15=(substr(titre,11,16))== "525-07" |  (substr(titre,11,16))== "525-08" |  (substr(titre,11,16))== "525-09" |  (substr(titre,11,16))== "525-10" |  (substr(titre,11,16))== "525-11" |  (substr(titre,11,16))== "525-12"
# x24=(substr(titre,11,16))== "525-13" |  (substr(titre,11,16))== "525-14" |  (substr(titre,11,16))== "525-15" |  (substr(titre,11,16))== "525-16" |  (substr(titre,11,16))== "525-17" |  (substr(titre,11,16))== "525-18"
# x7= (substr(titre,11,16))== "685-01" |  (substr(titre,11,16))== "685-02" |  (substr(titre,11,16))== "685-03" |  (substr(titre,11,16))== "685-04" |  (substr(titre,11,16))== "685-05" |  (substr(titre,11,16))== "685-06"
# x16=(substr(titre,11,16))== "685-07" |  (substr(titre,11,16))== "685-08" |  (substr(titre,11,16))== "685-09" |  (substr(titre,11,16))== "685-10" |  (substr(titre,11,16))== "685-11" |  (substr(titre,11,16))== "685-12"
# x25=(substr(titre,11,16))== "685-13" |  (substr(titre,11,16))== "685-14" |  (substr(titre,11,16))== "685-15" |  (substr(titre,11,16))== "685-16" |  (substr(titre,11,16))== "685-17" |  (substr(titre,11,16))== "685-18"
# x8= (substr(titre,11,16))== "747-01" |  (substr(titre,11,16))== "747-02" |  (substr(titre,11,16))== "747-03" |  (substr(titre,11,16))== "747-04" |  (substr(titre,11,16))== "747-05" |  (substr(titre,11,16))== "747-06"
# x17=(substr(titre,11,16))== "747-07" |  (substr(titre,11,16))== "747-08" |  (substr(titre,11,16))== "747-09" |  (substr(titre,11,16))== "747-10" |  (substr(titre,11,16))== "747-11" |  (substr(titre,11,16))== "747-12"
# x26=(substr(titre,11,16))== "747-13" |  (substr(titre,11,16))== "747-14" |  (substr(titre,11,16))== "747-15" |  (substr(titre,11,16))== "747-16" |  (substr(titre,11,16))== "747-17" |  (substr(titre,11,16))== "747-18"
# x9= (substr(titre,11,16))== "877-01" |  (substr(titre,11,16))== "877-02" |  (substr(titre,11,16))== "877-03" |  (substr(titre,11,16))== "877-04" |  (substr(titre,11,16))== "877-05" |  (substr(titre,11,16))== "877-06"
# x18=(substr(titre,11,16))== "877-07" |  (substr(titre,11,16))== "877-08" |  (substr(titre,11,16))== "877-09" |  (substr(titre,11,16))== "877-10" |  (substr(titre,11,16))== "877-11" |  (substr(titre,11,16))== "877-12"
# x27=(substr(titre,11,16))== "877-13" |  (substr(titre,11,16))== "877-14" |  (substr(titre,11,16))== "877-15" |  (substr(titre,11,16))== "877-16" |  (substr(titre,11,16))== "877-17" |  (substr(titre,11,16))== "877-18"

# #Comparaison de 222 et 787
# x1= (substr(titre,11,16))== "015-01" |  (substr(titre,11,16))== "015-02" |  (substr(titre,11,16))== "015-03" |  (substr(titre,11,16))== "015-04" |  (substr(titre,11,16))== "015-05" |  (substr(titre,11,16))== "015-06"
# x10=(substr(titre,11,16))== "015-07" |  (substr(titre,11,16))== "015-08" |  (substr(titre,11,16))== "015-09" |  (substr(titre,11,16))== "015-10" |  (substr(titre,11,16))== "015-11" |  (substr(titre,11,16))== "015-12"
# x19=(substr(titre,11,16))== "015-13" |  (substr(titre,11,16))== "015-14" |  (substr(titre,11,16))== "015-15" |  (substr(titre,11,16))== "015-16" |  (substr(titre,11,16))== "015-17" |  (substr(titre,11,16))== "015-18"
# x2= (substr(titre,11,16))== "169-01" |  (substr(titre,11,16))== "169-02" |  (substr(titre,11,16))== "169-03" |  (substr(titre,11,16))== "169-04" |  (substr(titre,11,16))== "169-05" |  (substr(titre,11,16))== "169-06"
# x11=(substr(titre,11,16))== "169-07" |  (substr(titre,11,16))== "169-08" |  (substr(titre,11,16))== "169-09" |  (substr(titre,11,16))== "169-10" |  (substr(titre,11,16))== "169-11" |  (substr(titre,11,16))== "169-12"
# x20=(substr(titre,11,16))== "169-13" |  (substr(titre,11,16))== "169-14" |  (substr(titre,11,16))== "169-15" |  (substr(titre,11,16))== "169-16" |  (substr(titre,11,16))== "169-17" |  (substr(titre,11,16))== "169-18"
# x3= (substr(titre,11,16))== "222-01" |  (substr(titre,11,16))== "222-02" |  (substr(titre,11,16))== "222-03" |  (substr(titre,11,16))== "222-04" |  (substr(titre,11,16))== "222-05" |  (substr(titre,11,16))== "222-06"
# x12=(substr(titre,11,16))== "222-07" |  (substr(titre,11,16))== "222-08" |  (substr(titre,11,16))== "222-09" |  (substr(titre,11,16))== "222-10" |  (substr(titre,11,16))== "222-11" |  (substr(titre,11,16))== "222-12"
# x21=(substr(titre,11,16))== "222-13" |  (substr(titre,11,16))== "222-14" |  (substr(titre,11,16))== "222-15" |  (substr(titre,11,16))== "222-16" |  (substr(titre,11,16))== "222-17" |  (substr(titre,11,16))== "222-18"
# x4= (substr(titre,11,16))== "471-01" |  (substr(titre,11,16))== "471-02" |  (substr(titre,11,16))== "471-03" |  (substr(titre,11,16))== "471-04" |  (substr(titre,11,16))== "471-05" |  (substr(titre,11,16))== "471-06"
# x13=(substr(titre,11,16))== "471-07" |  (substr(titre,11,16))== "471-08" |  (substr(titre,11,16))== "471-09" |  (substr(titre,11,16))== "471-10" |  (substr(titre,11,16))== "471-11" |  (substr(titre,11,16))== "471-12"
# x22=(substr(titre,11,16))== "471-13" |  (substr(titre,11,16))== "471-14" |  (substr(titre,11,16))== "471-15" |  (substr(titre,11,16))== "471-16" |  (substr(titre,11,16))== "471-17" |  (substr(titre,11,16))== "471-18"
# x5= (substr(titre,11,16))== "525-01" |  (substr(titre,11,16))== "525-02" |  (substr(titre,11,16))== "525-03" |  (substr(titre,11,16))== "525-04" |  (substr(titre,11,16))== "525-05" |  (substr(titre,11,16))== "525-06"
# x14=(substr(titre,11,16))== "525-07" |  (substr(titre,11,16))== "525-08" |  (substr(titre,11,16))== "525-09" |  (substr(titre,11,16))== "525-10" |  (substr(titre,11,16))== "525-11" |  (substr(titre,11,16))== "525-12"
# x23=(substr(titre,11,16))== "525-13" |  (substr(titre,11,16))== "525-14" |  (substr(titre,11,16))== "525-15" |  (substr(titre,11,16))== "525-16" |  (substr(titre,11,16))== "525-17" |  (substr(titre,11,16))== "525-18"
# x6= (substr(titre,11,16))== "685-01" |  (substr(titre,11,16))== "685-02" |  (substr(titre,11,16))== "685-03" |  (substr(titre,11,16))== "685-04" |  (substr(titre,11,16))== "685-05" |  (substr(titre,11,16))== "685-06"
# x15=(substr(titre,11,16))== "685-07" |  (substr(titre,11,16))== "685-08" |  (substr(titre,11,16))== "685-09" |  (substr(titre,11,16))== "685-10" |  (substr(titre,11,16))== "685-11" |  (substr(titre,11,16))== "685-12"
# x24=(substr(titre,11,16))== "685-13" |  (substr(titre,11,16))== "685-14" |  (substr(titre,11,16))== "685-15" |  (substr(titre,11,16))== "685-16" |  (substr(titre,11,16))== "685-17" |  (substr(titre,11,16))== "685-18"
# x7= (substr(titre,11,16))== "747-01" |  (substr(titre,11,16))== "747-02" |  (substr(titre,11,16))== "747-03" |  (substr(titre,11,16))== "747-04" |  (substr(titre,11,16))== "747-05" |  (substr(titre,11,16))== "747-06"
# x16=(substr(titre,11,16))== "747-07" |  (substr(titre,11,16))== "747-08" |  (substr(titre,11,16))== "747-09" |  (substr(titre,11,16))== "747-10" |  (substr(titre,11,16))== "747-11" |  (substr(titre,11,16))== "747-12"
# x25=(substr(titre,11,16))== "747-13" |  (substr(titre,11,16))== "747-14" |  (substr(titre,11,16))== "747-15" |  (substr(titre,11,16))== "747-16" |  (substr(titre,11,16))== "747-17" |  (substr(titre,11,16))== "747-18"
# x8= (substr(titre,11,16))== "787-01" |  (substr(titre,11,16))== "787-02" |  (substr(titre,11,16))== "787-03" |  (substr(titre,11,16))== "787-04" |  (substr(titre,11,16))== "787-05" |  (substr(titre,11,16))== "787-06"
# x17=(substr(titre,11,16))== "787-07" |  (substr(titre,11,16))== "787-08" |  (substr(titre,11,16))== "787-09" |  (substr(titre,11,16))== "787-10" |  (substr(titre,11,16))== "787-11" |  (substr(titre,11,16))== "787-12"
# x26=(substr(titre,11,16))== "787-13" |  (substr(titre,11,16))== "787-14" |  (substr(titre,11,16))== "787-15" |  (substr(titre,11,16))== "787-16" |  (substr(titre,11,16))== "787-17" |  (substr(titre,11,16))== "787-18"
# x9= (substr(titre,11,16))== "877-01" |  (substr(titre,11,16))== "877-02" |  (substr(titre,11,16))== "877-03" |  (substr(titre,11,16))== "877-04" |  (substr(titre,11,16))== "877-05" |  (substr(titre,11,16))== "877-06"
# x18=(substr(titre,11,16))== "877-07" |  (substr(titre,11,16))== "877-08" |  (substr(titre,11,16))== "877-09" |  (substr(titre,11,16))== "877-10" |  (substr(titre,11,16))== "877-11" |  (substr(titre,11,16))== "877-12"
# x27=(substr(titre,11,16))== "877-13" |  (substr(titre,11,16))== "877-14" |  (substr(titre,11,16))== "877-15" |  (substr(titre,11,16))== "877-16" |  (substr(titre,11,16))== "877-17" |  (substr(titre,11,16))== "877-18"

# #Comparaison 509 et 787
# x1= (substr(titre,11,16))== "015-01" |  (substr(titre,11,16))== "015-02" |  (substr(titre,11,16))== "015-03" |  (substr(titre,11,16))== "015-04" |  (substr(titre,11,16))== "015-05" |  (substr(titre,11,16))== "015-06"
# x10=(substr(titre,11,16))== "015-07" |  (substr(titre,11,16))== "015-08" |  (substr(titre,11,16))== "015-09" |  (substr(titre,11,16))== "015-10" |  (substr(titre,11,16))== "015-11" |  (substr(titre,11,16))== "015-12"
# x19=(substr(titre,11,16))== "015-13" |  (substr(titre,11,16))== "015-14" |  (substr(titre,11,16))== "015-15" |  (substr(titre,11,16))== "015-16" |  (substr(titre,11,16))== "015-17" |  (substr(titre,11,16))== "015-18"
# x2= (substr(titre,11,16))== "169-01" |  (substr(titre,11,16))== "169-02" |  (substr(titre,11,16))== "169-03" |  (substr(titre,11,16))== "169-04" |  (substr(titre,11,16))== "169-05" |  (substr(titre,11,16))== "169-06"
# x11=(substr(titre,11,16))== "169-07" |  (substr(titre,11,16))== "169-08" |  (substr(titre,11,16))== "169-09" |  (substr(titre,11,16))== "169-10" |  (substr(titre,11,16))== "169-11" |  (substr(titre,11,16))== "169-12"
# x20=(substr(titre,11,16))== "169-13" |  (substr(titre,11,16))== "169-14" |  (substr(titre,11,16))== "169-15" |  (substr(titre,11,16))== "169-16" |  (substr(titre,11,16))== "169-17" |  (substr(titre,11,16))== "169-18"
# x3= (substr(titre,11,16))== "471-01" |  (substr(titre,11,16))== "471-02" |  (substr(titre,11,16))== "471-03" |  (substr(titre,11,16))== "471-04" |  (substr(titre,11,16))== "471-05" |  (substr(titre,11,16))== "471-06"
# x12=(substr(titre,11,16))== "471-07" |  (substr(titre,11,16))== "471-08" |  (substr(titre,11,16))== "471-09" |  (substr(titre,11,16))== "471-10" |  (substr(titre,11,16))== "471-11" |  (substr(titre,11,16))== "471-12"
# x21=(substr(titre,11,16))== "471-13" |  (substr(titre,11,16))== "471-14" |  (substr(titre,11,16))== "471-15" |  (substr(titre,11,16))== "471-16" |  (substr(titre,11,16))== "471-17" |  (substr(titre,11,16))== "471-18"
# x4= (substr(titre,11,16))== "509-01" |  (substr(titre,11,16))== "509-02" |  (substr(titre,11,16))== "509-03" |  (substr(titre,11,16))== "509-04" |  (substr(titre,11,16))== "509-05" |  (substr(titre,11,16))== "509-06"
# x13=(substr(titre,11,16))== "509-07" |  (substr(titre,11,16))== "509-08" |  (substr(titre,11,16))== "509-09" |  (substr(titre,11,16))== "509-10" |  (substr(titre,11,16))== "509-11" |  (substr(titre,11,16))== "509-12"
# x22=(substr(titre,11,16))== "509-13" |  (substr(titre,11,16))== "509-14" |  (substr(titre,11,16))== "509-15" |  (substr(titre,11,16))== "509-16" |  (substr(titre,11,16))== "509-17" |  (substr(titre,11,16))== "509-18"
# x5= (substr(titre,11,16))== "525-01" |  (substr(titre,11,16))== "525-02" |  (substr(titre,11,16))== "525-03" |  (substr(titre,11,16))== "525-04" |  (substr(titre,11,16))== "525-05" |  (substr(titre,11,16))== "525-06"
# x14=(substr(titre,11,16))== "525-07" |  (substr(titre,11,16))== "525-08" |  (substr(titre,11,16))== "525-09" |  (substr(titre,11,16))== "525-10" |  (substr(titre,11,16))== "525-11" |  (substr(titre,11,16))== "525-12"
# x23=(substr(titre,11,16))== "525-13" |  (substr(titre,11,16))== "525-14" |  (substr(titre,11,16))== "525-15" |  (substr(titre,11,16))== "525-16" |  (substr(titre,11,16))== "525-17" |  (substr(titre,11,16))== "525-18"
# x6= (substr(titre,11,16))== "685-01" |  (substr(titre,11,16))== "685-02" |  (substr(titre,11,16))== "685-03" |  (substr(titre,11,16))== "685-04" |  (substr(titre,11,16))== "685-05" |  (substr(titre,11,16))== "685-06"
# x15=(substr(titre,11,16))== "685-07" |  (substr(titre,11,16))== "685-08" |  (substr(titre,11,16))== "685-09" |  (substr(titre,11,16))== "685-10" |  (substr(titre,11,16))== "685-11" |  (substr(titre,11,16))== "685-12"
# x24=(substr(titre,11,16))== "685-13" |  (substr(titre,11,16))== "685-14" |  (substr(titre,11,16))== "685-15" |  (substr(titre,11,16))== "685-16" |  (substr(titre,11,16))== "685-17" |  (substr(titre,11,16))== "685-18"
# x7= (substr(titre,11,16))== "747-01" |  (substr(titre,11,16))== "747-02" |  (substr(titre,11,16))== "747-03" |  (substr(titre,11,16))== "747-04" |  (substr(titre,11,16))== "747-05" |  (substr(titre,11,16))== "747-06"
# x16=(substr(titre,11,16))== "747-07" |  (substr(titre,11,16))== "747-08" |  (substr(titre,11,16))== "747-09" |  (substr(titre,11,16))== "747-10" |  (substr(titre,11,16))== "747-11" |  (substr(titre,11,16))== "747-12"
# x25=(substr(titre,11,16))== "747-13" |  (substr(titre,11,16))== "747-14" |  (substr(titre,11,16))== "747-15" |  (substr(titre,11,16))== "747-16" |  (substr(titre,11,16))== "747-17" |  (substr(titre,11,16))== "747-18"
# x8= (substr(titre,11,16))== "787-01" |  (substr(titre,11,16))== "787-02" |  (substr(titre,11,16))== "787-03" |  (substr(titre,11,16))== "787-04" |  (substr(titre,11,16))== "787-05" |  (substr(titre,11,16))== "787-06"
# x17=(substr(titre,11,16))== "787-07" |  (substr(titre,11,16))== "787-08" |  (substr(titre,11,16))== "787-09" |  (substr(titre,11,16))== "787-10" |  (substr(titre,11,16))== "787-11" |  (substr(titre,11,16))== "787-12"
# x26=(substr(titre,11,16))== "787-13" |  (substr(titre,11,16))== "787-14" |  (substr(titre,11,16))== "787-15" |  (substr(titre,11,16))== "787-16" |  (substr(titre,11,16))== "787-17" |  (substr(titre,11,16))== "787-18"
# x9= (substr(titre,11,16))== "877-01" |  (substr(titre,11,16))== "877-02" |  (substr(titre,11,16))== "877-03" |  (substr(titre,11,16))== "877-04" |  (substr(titre,11,16))== "877-05" |  (substr(titre,11,16))== "877-06"
# x18=(substr(titre,11,16))== "877-07" |  (substr(titre,11,16))== "877-08" |  (substr(titre,11,16))== "877-09" |  (substr(titre,11,16))== "877-10" |  (substr(titre,11,16))== "877-11" |  (substr(titre,11,16))== "877-12"
# x27=(substr(titre,11,16))== "877-13" |  (substr(titre,11,16))== "877-14" |  (substr(titre,11,16))== "877-15" |  (substr(titre,11,16))== "877-16" |  (substr(titre,11,16))== "877-17" |  (substr(titre,11,16))== "877-18"


sp1=sp[(x1==TRUE),] #1
sp2=sp[(x2==TRUE),] #2
sp3=sp[(x3==TRUE),] #3
sp4=sp[(x4==TRUE),] #4
sp5=sp[(x5==TRUE),] #5
sp6=sp[(x6==TRUE),] #6
sp7=sp[(x7==TRUE),] #7
sp8=sp[(x8==TRUE),] #8
sp9=sp[(x9==TRUE),] #9
sp10=sp[(x10==TRUE),] #10
sp11=sp[(x11==TRUE),] #11
sp12=sp[(x12==TRUE),] #12
sp13=sp[(x13==TRUE),] #13
sp14=sp[(x14==TRUE),] #14
sp15=sp[(x15==TRUE),] #15
sp16=sp[(x16==TRUE),] #16
sp17=sp[(x17==TRUE),] #17
sp18=sp[(x18==TRUE),] #18
sp19=sp[(x19==TRUE),] #19
sp20=sp[(x20==TRUE),] #20
sp21=sp[(x21==TRUE),] #21
sp22=sp[(x22==TRUE),] #22
sp23=sp[(x23==TRUE),] #23
sp24=sp[(x24==TRUE),] #24
sp25=sp[(x25==TRUE),] #25
sp26=sp[(x26==TRUE),] #26
sp27=sp[(x27==TRUE),] #27
sp28=sp[(x28==TRUE),] #28
sp29=sp[(x29==TRUE),] #29
sp30=sp[(x30==TRUE),] #30

# #Rajoute +30 a l'indice de la matrice pour la 2nd date
#   for (i in 1:30) {
#   commd=paste("sp",i+30,"=sp",i,sep="")
#   eval(parse(text=commd))
#   }


id=1:30
dim(id)=c(10,3)
id=t(id)

## Definition des matrices de resultat final
# Creation de la matrice des perok finale
perok_final=matrix(nrow = repet, ncol = 1)
perok_final_C=matrix(nrow = repet, ncol = 1)
perok_final_G=matrix(nrow = repet, ncol = 1)
perok_final_S=matrix(nrow = repet, ncol = 1)
perok_final_C_cep=matrix(nrow = repet, ncol = 1)
perok_final_G_cep=matrix(nrow = repet, ncol = 1)
perok_final_S_cep=matrix(nrow = repet, ncol = 1)
perok_final_cepages=matrix(nrow = repet, ncol = 1)
perok_finalm0=matrix(nrow = repet, ncol = ncmax)
perok_finalm0C=matrix(nrow = repet, ncol = ncmax)
perok_finalm0G=matrix(nrow = repet, ncol = ncmax)
perok_finalm0S=matrix(nrow = repet, ncol = ncmax)
## Creation matrice de % de mauvais classements par clone
mc=matrix(nrow = ncmax,ncol = c)

## Creation de la matrice des VL et perok maximaux
maxi_final=matrix(nrow= repet, ncol = 2)
maxi_finalC=matrix(nrow= repet, ncol = 2)
maxi_finalG=matrix(nrow= repet, ncol = 2)
maxi_finalS=matrix(nrow= repet, ncol = 2)
## Creation de la matrice de % de mauvais classements
mc_final=matrix(nrow= repet, ncol = length(levels(class)))
## Creation d'un matrice cubique pour enregistrer les tables de contingence
t_final=array(dim=c(c,c,repet))
## Noms des colonnes et des lignes
colnames(t_final)=c(basename(levels(class)))
rownames(t_final)=c(basename(levels(class)))
colnames(maxi_final)= c("maxi.id","perok max")
colnames(mc_final)= c(basename(levels(class)))

###s?paration validation calibration PLSDA###
#set.seed(1) # fixe le tirage aleatoire
for(j in 1:repet) {
#
# id2=matrix( ,nrow=1, ncol=10)
#  for (i in 1:3) {
#    id2[,i]=sample(id[,i],1, replace = FALSE)


 idval=matrix( ,nrow=1, ncol=10) #matrice validation PLSDA
 idcal=matrix( ,nrow=2, ncol=10) #matrice calibration PLSDA
 idcal2=matrix( ,nrow=2, ncol=10) #matrice calibration PLSDA
 # idval2=matrix( ,nrow=1, ncol=10) #idval pour CV
 # idcal2=matrix( ,nrow=1, ncol=10) #idcal pour CV

  for (i in 1:10) {
    icol=sample(1:3,1)
    # id2[,i]=sample(id[,i],1, replace = FALSE)
    idval[,i]=id[icol,i]
    idcal[,i]=id[c(1:3)[-icol],i]

    }

 for (i in 1:10) {
   # icol2=sample(1:2,1)
   idcal2[,i]=idcal[sample(1:2),i]
   # idval2[,i]=idcal[icol2,i]
   # idcal2[,i]=idcal[c(1:2)[-icol2],i]

     }

 ###PLSDA cepages

commd2=paste("sp_val=rbind(sp",idval[1,1],",sp",idval[1,2],",sp",idval[1,3],",sp",idval[1,4],",sp",idval[1,5],",sp",idval[1,6],",sp",idval[1,7],",sp",idval[1,8],",sp",idval[1,9],",sp",idval[1,10],")",sep="")
 #commd2=paste("sp_val=rbind(sp",idval[1,1],",sp",idval[1,2],",sp",idval[1,3],",sp",idval[1,4],",sp",idval[1,5],",sp",idval[1,6],",sp",idval[1,7],",sp",idval[1,8],",sp",idval[1,9],")",sep="") # Si on supprime 1 clone
 # commd2=paste("sp_val=rbind(sp",idval[1,1],",sp",idval[1,2],",sp",idval[1,3],",sp",idval[1,4],",sp",idval[1,5],",sp",idval[1,6],",sp",idval[1,7],",sp",idval[1,8],")",sep="")
 eval(parse(text=commd2))


 id_val=which(rownames(sp)  %in%  rownames(sp_val))

 ##On selectionne les spectres ayant ces num?ros dans le jeu de validation, les autres vont dans le jeu de calibration
 sp_val=sp[id_val,]
 sp_cal=sp[-id_val,]
 class_val=class[id_val]
 class_cal=class[-id_val]

 ## Creation des jeux d'apprentissage et validation
 predm0=as.data.frame(matrix(nrow = sum(iok2), ncol = ncmax))

## Boucle CV
   for (i in 1:k) {
    #commd=paste("sp_val=rbind(sp",idcal2[i,1],",sp",idcal2[i,2],",sp",idcal2[i,3],")",sep="")   # pour CS et G
    commd=paste("spvalCV=rbind(sp",idcal2[i,1],",sp",idcal2[i,2],",sp",idcal2[i,3],",sp",idcal2[i,4],")",sep="") # pour S
    #commd=paste("spvalCV=rbind(sp",idcal2[i,1],",sp",idcal2[i,2],",sp",idcal2[i,3],",sp",idcal2[i,4],",sp",idcal2[i,5],",sp",idcal2[i,6],",sp",idcal2[i,7],",sp",idcal2[i,8],",sp",idcal2[i,9],",sp",idcal2[i,10],")",sep="")
  # commd=paste("spvalCV=rbind(sp",idcal2[i,1],",sp",idcal2[i,2],",sp",idcal2[i,3],",sp",idcal2[i,4],",sp",idcal2[i,5],",sp",idcal2[i,6],",sp",idcal2[i,7],",sp",idcal2[i,8],")",sep="")
    # commd=paste("spvalCV=rbind(sp",idcal2[i,1],",sp",idcal2[i,2],",sp",idcal2[i,3],",sp",idcal2[i,4],",sp",idcal2[i,5],",sp",idcal2[i,6],",sp",idcal2[i,7],",sp",idcal2[i,8],",sp",idcal2[i,9],")",sep="")
     eval(parse(text=commd))

    idvalCV=which(rownames(sp_cal)  %in%  rownames(spvalCV))

    spvalCV=sp_cal[idvalCV,]       # matrice du jeu de validation
    class_valCV=class_cal[idvalCV]  #identifiants des classes du jeu de validation
    spcalCV=sp_cal[-idvalCV,]      #matrice du jeu de calibration compos?e de tout ce qui n'est pas en validation
    class_calCV=class_cal[-idvalCV] #identifiants des classes du jeu de calibration

    ## PLSDA and application to have loadings and scores
    rplsda=caret::plsda(spcalCV, class_calCV,ncomp=ncmax)
    sccalCV=rplsda$scores
    spvalCV_c=scale(spvalCV,center=rplsda$Xmeans,scale = F)
    scvalCV=spvalCV_c%*%rplsda$projection  # score_val=predict(rplsda,sc_val,type="scores") : ne marche pas

      for (ii in 2:ncmax) {

     ## Validation
     predm0[idvalCV,ii]=SIGNE_maha0(sccalCV[,1:ii], class_calCV, scvalCV[,1:ii])$class
       }
   }

 ## Table de contingence CV
  tsm0=lapply(as.list(predm0), class, FUN = table)

 ## Matrice mauvais classements par clone CV
  diagsm0=lapply(tsm0, FUN = diag)

 ## Pourcentage de bien classes CV
  perokm0=100*unlist(lapply(diagsm0, FUN = sum))/length(class)

 ## Pourcentage de bien classes CV
  maxi=max(perokm0)
  maxi.id=which.max(perokm0)

 ### Enregistrement des matrices de resultat final
 ##Remplissage de la matrice des perok finale
 perok_finalm0[j,]=perokm0

 ## Remplissage de la VL max et de son % de bons classements globaux
 maxi_final[j,1]=maxi.id
 maxi_final[j,2]=maxi
 ## Remplissage de la matrice de mauvais classements par clone
 mc_final[j,]=mc[maxi.id,]


###PLSDA on Maha scores
## Calibration
rplsda=caret::plsda(sp_cal, class_cal,ncomp=10)
sc_cal=rplsda$scores

## Validation
sp_val_c=scale(sp_val,center=rplsda$Xmeans,scale = F)
sc_val=sp_val_c%*%rplsda$projection
res_val=SIGNE_maha0(sc_cal[,1:10], class_cal, sc_val[,1:10])$class

cepage=table (res_val,class_val)
#print (cepage)


###En fonction des clones

## Creation de la matrice de classes clones
class_clones=as.factor(substr(rownames(sp),11,13))
## Variable qui mesure le nombre de classes
c=length(levels(class_clones))

## Separation de sp_cal en 3 jeux de calibration par cepage
aC= substr(rownames(sp_cal),18,18)=="C"
aG= substr(rownames(sp_cal),18,18)=="G"
aS= substr(rownames(sp_cal),18,18)=="S"

sp_cal_C =sp_cal[(aC==TRUE),]
sp_cal_G =sp_cal[(aG==TRUE),]
sp_cal_S =sp_cal[(aS==TRUE),]

###Identifiants des matrices de calibration
##Cabernet
id_cal_C=which(rownames(sp)  %in%  rownames(sp_cal_C))
class_cal_C=droplevels(class_clones[id_cal_C])

##Gamay
id_cal_G=which(rownames(sp)  %in%  rownames(sp_cal_G))
class_cal_G=droplevels(class_clones[id_cal_G])

##Syrah
id_cal_S=which(rownames(sp)  %in%  rownames(sp_cal_S))
class_cal_S=droplevels(class_clones[id_cal_S])

# ###CV sur clones
# ##Formation de idcalC (Cabernet)
#   c1=idcal[1,1]
#   c2=idcal[1,2]
#   c3=idcal[1,3]
#   c4=idcal[2,1]
#   c5=idcal[2,2]
#   c6=idcal[2,3]
#
# idcalC=c(c1,c2,c3,c4,c5,c6)
# dim(idcalC)=c(3,2)
# idcalC=t(idcalC)
#
# idcal2C=matrix( ,nrow=2, ncol=3)
#
# for (i in 1:3) {
#   idcal2C[,i]=idcalC[sample(1:2),i]
# }
#
# ## Boucle CV (Cabernet)
# for (i in 1:k) {
#   commdC=paste("spvalCVC=rbind(sp",idcal2C[i,1],",sp",idcal2C[i,2],",sp",idcal2C[i,3],")",sep="")
#   eval(parse(text=commdC))
#
#   idvalCVC=which(rownames(sp_cal_C)  %in%  rownames(spvalCVC))
#
#   spvalCVC=sp_cal_C[idvalCVC,]       # matrice du jeu de validation
#   class_valCVC=droplevels(class_cal_C[idvalCVC])  #identifiants des classes du jeu de validation
#   spcalCVC=sp_cal_C[-idvalCVC,]      #matrice du jeu de calibration compos?e de tout ce qui n'est pas en validation
#   class_calCVC=droplevels(class_cal_C[-idvalCVC]) #identifiants des classes du jeu de calibration
#
#   ## PLSDA and application to have loadings and scores (Cabernet)
#   rplsdaCVC=caret::plsda(spcalCVC, class_calCVC,ncomp=ncmax)
#   sccalCVC=rplsdaCVC$scores
#   spvalCVC_c=scale(spvalCVC,center=rplsdaCVC$Xmeans,scale = F)
#   scvalCVC=spvalCVC_c%*%rplsdaCVC$projection  # score_val=predict(rplsda,sc_val,type="scores") : ne marche pas
#
#   ## Creation des jeux d'apprentissage et validation
#   iok3=substr(rownames(sp_cal_C),1,9) %in% dates
#   predm0C=as.data.frame(matrix(nrow = sum(iok3), ncol = ncmax))
#
#
#   for (ii in 2:ncmax) {
#
#     ## Validation
#     predm0C[idvalCVC,ii]=SIGNE_maha0(sccalCVC[,1:ii], class_calCVC, scvalCVC[,1:ii])$class
#   }
# }
#
#
# ## Table de contingence CV
# tsm0C=lapply(as.list(predm0C), class_cal_C, FUN = table)
#
#
# ## Matrice mauvais classements par clone CV
# diagsm0C=lapply(tsm0C, FUN = diag)
#
# ## Pourcentage de bien classes CV
# perokm0C=100*unlist(lapply(diagsm0C, FUN = sum))/length(class_valCVC)
#
# ## Pourcentage de bien classes CV
# maxiC=max(perokm0C)
# maxi.idC=which.max(perokm0C)
#
# ### Enregistrement des matrices de resultat final
# ##Remplissage de la matrice des perok finale
# perok_finalm0C[j,]=perokm0C
#
# ## Remplissage de la VL max et de son % de bons classements globaux
# maxi_finalC[j,1]=maxi.idC
# maxi_finalC[j,2]=maxiC
#
# ##Formation de idcalG (Gamay)
#   g1=idcal[1,4]
#   g2=idcal[1,5]
#   g3=idcal[1,6]
#   g4=idcal[2,4]
#   g5=idcal[2,5]
#   g6=idcal[2,6]
#
# idcalG=c(g1,g2,g3,g4,g5,g6)
# dim(idcalG)=c(3,2)
# idcalG=t(idcalG)
#
# idcal2G=matrix( ,nrow=2, ncol=3)
#
# for (i in 1:3) {
#   idcal2G[,i]=idcalG[sample(1:2),i]
# }
#
# ## Boucle CV (Gamay)
# for (i in 1:k) {
#   commdG=paste("spvalCVG=rbind(sp",idcal2G[i,1],",sp",idcal2G[i,2],",sp",idcal2G[i,3],")",sep="")
#   eval(parse(text=commdG))
#
#   idvalCVG=which(rownames(sp_cal_G)  %in%  rownames(spvalCVG))
#
#   spvalCVG=sp_cal_G[idvalCVG,]       # matrice du jeu de validation
#   class_valCVG=droplevels(class_cal_G[idvalCVG])  #identifiants des classes du jeu de validation
#   spcalCVG=sp_cal_G[-idvalCVG,]      #matrice du jeu de calibration compos?e de tout ce qui n'est pas en validation
#   class_calCVG=droplevels(class_cal_G[-idvalCVG]) #identifiants des classes du jeu de calibration
#
#   ## PLSDA and application to have loadings and scores (Gamay)
#   rplsdaCVG=caret::plsda(spcalCVG, class_calCVG,ncomp=ncmax)
#   sccalCVG=rplsdaCVG$scores
#   spvalCVG_c=scale(spvalCVG,center=rplsdaCVG$Xmeans,scale = F)
#   scvalCVG=spvalCVG_c%*%rplsdaCVG$projection  # score_val=predict(rplsda,sc_val,type="scores") : ne marche pas
#
#   ## Creation des jeux d'apprentissage et validation
#   iok4=substr(rownames(sp_cal_G),1,9) %in% dates
#   predm0G=as.data.frame(matrix(nrow = sum(iok4), ncol = ncmax))
#
#   for (ii in 2:ncmax) {
#
#     ## Validation
#     predm0G[idvalCVG,ii]=SIGNE_maha0(sccalCVG[,1:ii], class_calCVG, scvalCVG[,1:ii])$class
#   }
# }
#
# ## Table de contingence CV
# tsm0G=lapply(as.list(predm0G), class_cal_G, FUN = table)
#
# ## Matrice mauvais classements par clone CV
# diagsm0G=lapply(tsm0G, FUN = diag)
#
# ## Pourcentage de bien classes CV
# perokm0G=100*unlist(lapply(diagsm0G, FUN = sum))/length(class_valCVG)
#
# ## Pourcentage de bien classes CV
# maxiG=max(perokm0G)
# maxi.idG=which.max(perokm0G)
#
# ### Enregistrement des matrices de resultat final
# ##Remplissage de la matrice des perok finale
# perok_finalm0G[j,]=perokm0G
#
# ## Remplissage de la VL max et de son % de bons classements globaux
# maxi_finalG[j,1]=maxi.idG
# maxi_finalG[j,2]=maxiG

# ##Formation de idcalS (Syrah)
#   s1=idcal[1,7]
#   s2=idcal[1,8]
# #  s3=idcal[1,9]
# #  s4=idcal[1,10]
#   s5=idcal[2,7]
#   s6=idcal[2,8]
# #  s7=idcal[2,9]
# #  s8=idcal[2,10]
#
# idcalS= c(s1,s2,s5,s6)#c(s1,s2,s3,s4,s5,s6,s7,s8)
# dim(idcalS)=c(2,2)
# idcalS=t(idcalS)
#
# idcal2S=matrix( ,nrow=2, ncol=2)
#
# for (i in 1:2) {
#   idcal2S[,i]=idcalS[sample(1:2),i]
# }

# ## Boucle CV  (Syrah)
# for (i in 1:k) {
#   #commdS=paste("spvalCVS=rbind(sp",idcal2S[i,1],",sp",idcal2S[i,2],",sp",idcal2S[i,3],",sp",idcal2S[i,4],")",sep="")
#   commdS=paste("spvalCVS=rbind(sp",idcal2S[i,1],",sp",idcal2S[i,2],")",sep="")
#   eval(parse(text=commdS))
#
#   idvalCVS=which(rownames(sp_cal_S)  %in%  rownames(spvalCVS))
#
#   spvalCVS=sp_cal_S[idvalCVS,]       # matrice du jeu de validation
#   class_valCVS=droplevels(class_cal_S[idvalCVS])  #identifiants des classes du jeu de validation
#   spcalCVS=sp_cal_S[-idvalCVS,]      #matrice du jeu de calibration compos?e de tout ce qui n'est pas en validation
#   class_calCVS=droplevels(class_cal_S[-idvalCVS]) #identifiants des classes du jeu de calibration
#
#   ## PLSDA and application to have loadings and scores (Syrah)
#   rplsdaCVS=caret::plsda(spcalCVS, class_calCVS,ncomp=ncmax)
#   sccalCVS=rplsdaCVS$scores
#   spvalCVS_c=scale(spvalCVS,center=rplsdaCVS$Xmeans,scale = F)
#   scvalCVS=spvalCVS_c%*%rplsdaCVS$projection  # score_val=predict(rplsda,sc_val,type="scores") : ne marche pas
#
#   ## Creation des jeux d'apprentissage et validation
#    iok5=substr(rownames(sp_cal_S),1,9) %in% dates
#    predm0S=as.data.frame(matrix(nrow = sum(iok5), ncol = ncmax))
#
#   for (ii in 2:ncmax) {
#
#     ## Validation
#     predm0S[idvalCVS,ii]=SIGNE_maha0(sccalCVS[,1:ii], class_calCVS, scvalCVS[,1:ii])$class
#   }
# }

# ## Table de contingence CV
# tsm0S=lapply(as.list(predm0S), class_cal_S, FUN = table)
#
# ## Matrice mauvais classements par clone CV
# diagsm0S=lapply(tsm0S, FUN = diag)
#
# ## Pourcentage de bien classes CV
# perokm0S=100*unlist(lapply(diagsm0S, FUN = sum))/length(class_valCVS)
#
# ## Pourcentage de bien classes CV
# maxiS=max(perokm0S)
# maxi.idS=which.max(perokm0S)
#
# ### Enregistrement des matrices de resultat final
# ##Remplissage de la matrice des perok finale
# perok_finalm0S[j,]=perokm0S
#
# ## Remplissage de la VL max et de son % de bons classements globaux
# maxi_finalS[j,1]=maxi.idS
# maxi_finalS[j,2]=maxiS

###PLSDA sur clones sans CV
## Separation de sp_val en 3 jeux de validation par cepage
## Seulement avec les biens classes en cepages
test=cbind(res_val,class_val)
rownames(test)=rownames(sp[id_val,])

z1=(substr(rownames(test),18,18))=="C"
z2=(substr(rownames(test), 18, 18))=="G"
z3=(substr(rownames(test),18,18))=="S"

test1 =test[(z1==TRUE),]
test2 =test[(z2==TRUE),]
test3 =test[(z3==TRUE),]

test1a = test1[test1[,1]==1,]
test2a = test2[test2[,1]==2,]
test3a = test3[test3[,1]==3,]

sp_val_C =sp_val[rownames(test1a),]
sp_val_G =sp_val[rownames(test2a),]
sp_val_S =sp_val[rownames(test3a),]

##Avec tous
# bC= substr(rownames(sp_val),18,18)=="C"
# bG= substr(rownames(sp_val),18,18)=="G"
# bS= substr(rownames(sp_val),18,18)=="S"
#
# sp_val_C =sp_val[(bC==TRUE),]
# sp_val_G =sp_val[(bG==TRUE),]
# sp_val_S =sp_val[(bS==TRUE),]

###Selection des spectres
##Cabernet Sauvignon
id_val_C=which(rownames(sp)  %in%  rownames(sp_val_C))
class_val_C=droplevels(class_clones[id_val_C])

##Gamay
id_val_G=which(rownames(sp)  %in%  rownames(sp_val_G))
class_val_G=droplevels(class_clones[id_val_G])


##Syrah
id_val_S=which(rownames(sp)  %in%  rownames(sp_val_S))
class_val_S=droplevels(class_clones[id_val_S])

####PLSDA on Maha scores
### Calibration
## Cabernet Sauvignon
rplsdaC=caret::plsda(sp_cal_C, class_cal_C,ncomp= 8)# Modifier ncmax en fonction des resultats de CV_2
sc_cal_C=rplsdaC$scores

##Gamay
rplsdaG=caret::plsda(sp_cal_G, class_cal_G,ncomp=8)# Modifier ncmax en fonction des resultats de CV_2
sc_cal_G=rplsdaG$scores

##Syrah
rplsdaS=caret::plsda(sp_cal_S, class_cal_S,ncomp=8)# Modifier ncmax en fonction des resultats de CV_2
sc_cal_S=rplsdaS$scores

### Validation
## Cabernet Sauvignon
sp_val_c_C=scale(sp_val_C,center=rplsdaC$Xmeans,scale = F)
sc_val_C=sp_val_c_C%*%rplsdaC$projection
res_val_C=SIGNE_maha0(sc_cal_C[,1:8], class_cal_C, sc_val_C[,1:8])$class

Cabernet=table (res_val_C,class_val_C)
#print (Cabernet)

##Gamay
sp_val_c_G=scale(sp_val_G,center=rplsdaG$Xmeans,scale = F)
sc_val_G=sp_val_c_G%*%rplsdaG$projection
res_val_G=SIGNE_maha0(sc_cal_G[,1:8], class_cal_G, sc_val_G[,1:8])$class

Gamay=table (res_val_G,class_val_G)
#print (Gamay)

##Syrah
sp_val_c_S=scale(sp_val_S,center=rplsdaS$Xmeans,scale = F)
sc_val_S=sp_val_c_S%*%rplsdaS$projection
res_val_S=SIGNE_maha0(sc_cal_S[,1:8], class_cal_S, sc_val_S[,1:8])$class

Syrah=table (res_val_S,class_val_S)
#print (Syrah)

###Calcul des pourcentages de bons classements en fonction du tirage
##Somme des clones biens classes
clones_bc= sum(diag(Cabernet))+sum(diag(Gamay))+sum(diag(Syrah))
#print(clones_bc)

#Somme colonne table contingence cepages
sumcolcep=apply(cepage,2,sum)

##Nombre total de clones
total=clones_bc+(sum(cepage)-sum(diag(cepage)))+(sum(Cabernet)-sum(diag(Cabernet)))+(sum(Gamay)-sum(diag(Gamay)))+(sum(Syrah)-sum(diag(Syrah)))
total_C=sum(Cabernet)
total_G=sum(Gamay)
total_S=sum(Syrah)

total_C_cep=sum(diag(Cabernet))+(sum(Cabernet)-sum(diag(Cabernet)))+(sumcolcep[1]-cepage[1,1])
total_G_cep=sum(diag(Gamay))+(sum(Gamay)-sum(diag(Gamay)))+(sumcolcep[2]-cepage[2,2])
total_S_cep=sum(diag(Syrah))+(sum(Syrah)-sum(diag(Syrah)))+(sumcolcep[3]-cepage[3,3])

#print(total)

## Pourcentage total de clones biens classes(prise en compte erreur cepages)
perok=100*(clones_bc/total)
#print(perok)
perok_final[j,]=perok

## Pourcentage de clones biens classes
perok_C=100*(sum(diag(Cabernet))/total_C)
perok_G=100*(sum(diag(Gamay))/total_G)
perok_S=100*(sum(diag(Syrah))/total_S)

perok_final_C[j,]=perok_C
perok_final_G[j,]=perok_G
perok_final_S[j,]=perok_S



## Pourcentage de clones biens classes (prise en compte erreur cepages)
perok_C_cep=100*(sum(diag(Cabernet))/total_C_cep)
perok_G_cep=100*(sum(diag(Gamay))/total_G_cep)
perok_S_cep=100*(sum(diag(Syrah))/total_S_cep)

perok_final_C_cep[j,]=perok_C_cep
perok_final_G_cep[j,]=perok_G_cep
perok_final_S_cep[j,]=perok_S_cep


##Pourcentage de cepages biens classes
perok_cepages=100*(sum(diag(cepage))/sum(cepage))
#print(perok_cepages)
perok_final_cepages[j,]=perok_cepages

}

# print(perok_final)
# print(perok_final_cepages)
print(mean(perok_final))
print(mean(perok_final_cepages))
print(mean(perok_final_C))
print(mean(perok_final_G))
print(mean(perok_final_S))
print(mean(perok_final_C_cep))
print(mean(perok_final_G_cep))
print(mean(perok_final_S_cep))

###Sorties graphiques
## Tracage de l'evolution des perok en fonction du nombre de VL utilisees (cepages)
plot(colMeans(perok_finalm0), xlab= "Nombre de VL", ylab = "Pourcentage de biens class?s",pch=19, cex=1.5)
legend(ncmax*2/3,15,legend=c("Maha on PLSDA scores"),
       col=c("black"), lty=1, cex=0.8)

# ## Tracage de l'evolution des perok en fonction du nombre de VL utilisees (Cabernet)
# plot(colMeans(perok_finalm0C), xlab= "Nombre de VL", ylab = "Pourcentage de biens class?s",pch=19, cex=1.5)
# legend(ncmax*2/3,15,legend=c("Maha on PLSDA scores"),
#        col=c("black"), lty=1, cex=0.8)
#
# ## Tracage de l'evolution des perok en fonction du nombre de VL utilisees (Gamay)
# plot(colMeans(perok_finalm0G), xlab= "Nombre de VL", ylab = "Pourcentage de biens class?s",pch=19, cex=1.5)
# legend(ncmax*2/3,15,legend=c("Maha on PLSDA scores"),
#        col=c("black"), lty=1, cex=0.8)
#
# ## Tracage de l'evolution des perok en fonction du nombre de VL utilisees (Syrah)
# plot(colMeans(perok_finalm0S), xlab= "Nombre de VL", ylab = "Pourcentage de biens class?s",pch=19, cex=1.5)
# legend(ncmax*2/3,15,legend=c("Maha on PLSDA scores"),
#        col=c("black"), lty=1, cex=0.8)

## Clones
plot(perok_final, type="o",xlab= "Nombre de tirages", ylab = "Pourcentage de clones biens class?s",pch=19, cex=2)

##Cepages
plot(perok_final_cepages, type="o",xlab= "Nombre de tirages", ylab = "Pourcentage de c?pages biens class?s",pch=21, cex=2,bg="blue")
martinEcarnot/signe documentation built on Nov. 17, 2022, 1:49 a.m.

R Package Documentation

rdrr.io home R language documentation Run R code online

Browse R Packages

CRAN packages Bioconductor packages R-Forge packages GitHub packages

We want your feedback!

Note that we can't provide technical support on individual packages. You should contact the package authors for that.
GitHub issue tracker
ian@mutexlabs.com
Personal blog

 
Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.

Close